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‘ 'Atmospheric  density  profiles  are  presented  for  the  following  species; 
atomic  oxygen,  molecular  oxygen  in  theU'g  state,  ozone,  nitric  oxide, 
atomic  nitrogen,  nitrogen  dioxide,  carbon  dioxide,  and  water  vapor.  The 
profiles  are  intended  to  represent  mean  densities  at  midlatitudes  for  noon 
and  midnight  conditions.  A review  and  limited  discussion  of  recent  obser- 
vations and  calculations  are  given. 
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18.  SUPPLEMENTARY  NOTES  (Continued) 

EDITORS'  NOTE 


Volumes  13  to  17  were  originally  published  by  SAI  to  describe 
the  atmospheric,  geomagnetic,  and  high-altitude  energy  deposition  and 
neutral  heave  models  for  ROSCOE.  This  whole  section  of  code,  when 
associated  with  an  appropriate  DRIVER  subroutine,  operated  as  a package 
that  ran  independently  of  the  rest  of  the  ROSCOE  structure.  Provision 
was  also  made,  within  this  high-altitude  package,  for  two  completely 
independent  descriptions  of  atmospheric  heave,  each  with  its  own 
description  of  atmospheric  chemistry. 

When  GRC  incorporated  this  section  of  code  within  the  ROSCOE  i 
framework,  some  mcxlifications  were  necessary,  which  means  that  j 

some  of  the  descriptions  in  Volumes  13  to  17  are  inappropriate  to  ROSCOE  ^ 
as  it  now  exists.  In  particular,  the  NRL  heave  routines  (deck  NRLHYD)  i 
and  associated  chemistry  (deck  NRLCHM)  are  not  presently  used  in  j 

ROSCOE.  Three  other  subroutines  are  different;  subroutines  ATMOSU. 

EIF,  and  XTCOEF  correspond  to  the  ROSCOE  subroutines  ATMOS,  | 

EXPINT,  and  WDXP  respectively.  With  these  exceptions,  the  subroutines  | 
described  in  Volumes  13  to  17  correspond  exactly  to  those  currently  in  ^ 

ROSCOE. 

I 

I 

i 

The  bulk  of  this  work  ( 85  percent')  was  supported  by  the  author;  I 
the  remaining  portion  was  sponsored  by  the  Defense  Nuclear  Agency  (see 
Item  8 ) . 
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SUMMARY 


Atmosphei’ic  density  profiles  are  presented  for  the  following 
species:  atomic  oxygen,  molecular  oxygen  in  the  state,  ozone, 

nitric  oxide,  atomic  nitrogen,  nitrogen  dioxide,  carbon  dioxide  and 
water  vapor.  These  profiles  are  intended  to  represent  mean  densities 
at  midlatitudes  for  noon  and  midnight  conditions.  A review  and  limited 
discussion  of  recent  observations  and  calculations  are  given. 
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1.  INTRODUCTION 


Profiles  of  selected  atmospheric  species  have  been  constmcted 
based  on  a collation  of  observed  and  calculated  densities  from  ground 
level  to  several  hundred  kilometers  altitude.  The  profiles  are  intended 
to  represent  mean  densities  at  midlatitudes  primarily  in  the  noidhern 
hemisphere.  The  profiles  are  presented  for  only  two  diurnal  conditions: 
noon  and  midnight.  In  the  data  from  which  the  profiles  were  drawn, 
there  are  to  be  found  significant  variations  as  a result  of  geographical, 
seasonal,  solar  and  other  environmental  factors  but  the  intended  appb' 
cations  of  the  profiles  were  envisaged  to  be  not  so  precise  as  to  require 
the  incorporation  of  such  variations,  were  this  possible  with  the  present 
state  of  knowledge  about  each  species.  Furthermore,  there  is  consider- 
able scatter  in  obsei’vational  data  (as  well  as  in  calculated  densities)  in 
certain  ranges  of  altitude  which  may  result  not  only  from  exjieri mental 
difficulties  but  also  from  a lack  of  understanding  or  analysis  of  the  effect 
of  the  factors  mentioned  above. 

Onlv  mean  profiles  have  been  presenfed;  the  uncertainty  asso- 
ciated with  these  profiles  has  not  lieen  explicitly  given.  However,  a 
general  appreciation  of  the  uncertainties  in  the  mean  profiles  can  be  ob- 
tained from  the  collations  of  the  data  reviewed  in  this  report. 

Profiles  are  presented  for  the  following  species:  atomic  oxygen, 
molecular  oxygen  in  the  a^A^^,  state,  ozone,  nitric  oxide,  atomic  nit rogen, 
nitrogen  dioxide,  carbon  dioxide,  and  water  vapor.  This  list  is  re- 
stricted, of  course,  and  in  certain  altitude  ranges,  excepting  molecular 
oxygen  and  nitrogen,  important  species  are  not  included  as  is  the  case 
for  nitric  acid  in  the  troposphere  and  stratosphere. 


PAGE  MOT  FixMEi) 

HLAliiK 


Following  the  completion  of  the  profile  selection,  additional, 
important  data  have  appeared  in  the  literature.  These  data  are  presented 
separately  as  late  reports  in  the  sections  for  the  corresponding  species. 
Also,  note  that  the  survey  of  the  literature  was  not  intended  to  be  ex- 
haustive but  to  include  primarily  the  most  recent  obseiwations  and 
calculations. 
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2.  DENSITY  PROFILE  OF  ATOMIC  OXYGEN 


2.  1 INTRODUCTION 

An  atomic  oxytJ,en  profile  is  selected  for  altitudes  from  ground 
level  to  120  km  and  for  midlatitudes.  The  selection  is  based  on  obser- 
vations and  model  calculations.  At  present,  the  reported  obsen'ations 
of  atomic  oxvgen  are  restricted  to  altitudes  above  90  km;  l^elow  90  km, 
reliance  has  to  be  placed  on  model  calculations. 

Several  difficulties  are  encountered  when  attempting  to  select 
an  atomic  oxygen  profile.  Foi  altitudes  above  90  km,  the  obseiwations 
demonst rate  diu rnal  jDa-72d,  OD-73a,  Of-74,  Po-72a|,  seasonal  [DG-73, 
OD-73a,  Of-74,  Sc-74a,  SL-72aj,  and  latitudinal  and  longitudinal  [Da-72d. 
DG-73,  DG-74]  variations  in  the  atomic  oxygen  densitv:  these  variations 
cannot  yet  be  taken  into  account  in  a consistent  manner  as  there  are  im- 
portant problems  concerning  t ransport  of  atomic  o.xygen  which  have  vet 
to  be  understood.  In  addition,  the  obseiwations  are  sporadic,  thus  far, 
and  this  contributes  to  uncertainties  in  the  selected  profiles.  Below  90 
km,  where  model  calculations  are  the  source  of  atomic  oxvgen  densities, 
problems  involving  transport  coefficients  (i.  e.  , eddy  diffusion  coeffi- 
cients), boundary  conditions,  restricted  modelling  of  mass  motions  and 
simplified  sets  of  reactions  remain.  Adequate  comparisons  between 
different  model  calculations  have  not  been  made. 

2.  2 ATOMIC  OXYGEN  DENSITIES  ABOVE  90  KM 
2.  2.  1 Measu  re  merits 

The  atomic  o.xvgen  densitie.s  above  about  90  km  are  plotted  in 
Fig.  2-1.  Measurements  [NS-73,  OG-73.  Po-72a  ] for  altitudes  above 
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120  km  are  ^iven  on  Fig.  2-1  l)ut  will  not  be  discussed  in  general.  Con- 
sider first,  the  observational  data  for  rocket-based  experiments. 
Henderson  lHe-74a|  has  measured  atomic  oxygen  densities  by  use  of  a 
thin,  silver  film  whose  resistance  changes  in  the  event  that  oxygen 
atoms  are  adsorbed:  his  measurements  apply  to  a latitude  of  39  N. 

Both  the  original  and  corrected  [He-74a]  data  of  Henderson  are  plotted 
in  Fig.  2-1;  data  between  90  and  105  km  were  reduced  in  value  by 
Henderson  so  as  to  obtain  a reasonable,  composite  profile  over  the  alti- 
tude range  90  to  125  km.  The  data  of  Henderson  peak  near  97  km  and 
decline  rapidly  at  other  altitudes.  The  error  associated  with  Henderson's 
measurements  is  about  ±50%.  Measurements  were  also  made  |He-74a] 
at  a latitude  of  71  °N  but  these  measurements  are  not  included  in  Fig.  2-1 
as  our  present  interests  are  restricted  to  midlatitudes.  However,  note 
that  at  a latitude  of  71 ’N,  the  measured  oxygen  atom  densities  are  an 
order  of  magnitude  smaller  than  at  39  N,  using  the  same  experimental 
technique  for  the  range  of  altitudes  between  90  and  115  km,  and  become 
similar  only  above  120  km. 

An  extensive  atomic  oxygen  profile  has  been  obtained  by 
Offermann  et  al.  [OD-73a,  OG-73  | between  90  and  240  km  by  combining 
photometer  measurements  at  5577  A for  altitudes  below  118  km  with 
mass  spectrometer  measurements  at  altitudes  above  118  km.  The  mass 
spectrometer  walls  nere  cryogenically  cooled  to  8 K to  minimize  losses 
of  atomic  o.xygen  in  the  ion  source.  The  measurements  were  consistent 
in  that  the  photometer  and  mass  spectrometer  yielded  atomic  oxygen 
densities  which  were  in  good  agreement  in  the  overlap  altitude  range 
near  118  km.  The  aensity  of  atomic  oxvgen  peaks  at  about  97  km  and 
is  in  agreement  with  the  corrected  data  of  Henderson  |He-74a|  l^etween 
90  and  115  km.  Note  also,  the  good  agreement  above  120  km  between 
the  atomic  o.xvgen  densities  of  Offermann  and  Grossm.ann  [OG-73  | and 


those  of  the  CIRA  72  model  atmosphere.  Offermann  and  Drescher 
(OD-73aJ  summarize  previous  atomic  oxygen  measurements  (based 
mainly  on  photometer  results),  and  find  valuations  with  day  of  the  year, 
local  time  and  a weak  variation  with  magnetic  index. 

Scholz  and  Offermann  |SO-74b]  have  extended  the  mass  spec- 
trometer measurements  to  the  altitude  range  85  to  115  km  at  latitude 
39.6  N;  these  results  are  shown  in  Fig.  2-1;  their  oxygen  atom  densi- 
ties are  relatively  high,  particularly  at  altitudes  above  105  km.  The 
results  of  Offermann  and  colleagues  have  been  compared  with  previous 
mass  spectrometer  measurements  [SO-74b,  Of-74]  and  the  conclusion 
is  reached  that,  with  a cei’tain  degree  of  confidence,  the  previous  results 
obtained  with  uncooled  mass  spectrometer  walls  should  be  increased  by 
a factor  of  2.  5.  A review  of  the  older  mass  spectrometer  measure- 
ments has  not  been  included  in  the  present  review. 

Also  shown  in  Fig.  2-1,  is  the  range  of  atomic  oxygen  densi- 
ties at  95  km  deduced  from  photometer  measurements  of  the  5577  A air- 
glow  [OD-73a]. 

Measurements  of  the  atomic  oxv'gen  density  have  been  made  by 
Golomb  and  Good  |GG-72]  based  on  the  O-NO  chemiluminescence. 

These  measurements,  which  are  shown  on  Fig.  2-1,  were  made  at 
22  N and  59  N;  the  error  of  these  measurements  is  a factor  of  : 2.  5 
below  120  km.  For  the  measurements  at  59  N,  data  were  obtained  at 
224  5 hrs  (and  0300  hrs  but  these  measurements  are  not  represented  on 
Fig.  2-1).  Golomb  and  Good  note  that  the  peak  of  the  oxygen  atom  den- 
sity seems  to  ascend  from  about  93  to  99  km  as  the  night  progresses  and 
this  can  be  understood  from  their  measurements  and  also  bv  comparison 
with  Henderson's  (at  0230  hr)  if  considerations  of  latitude  are  neglected. 
The  results  of  Golomb  and  Good  between  110  and  150  km  are  not  shown: 


14 


in  this  altitude  region,  they  found  a minimum  at  120  km  in  the  atomic 

oxygen  profiles  with  values  of  the  density  ranging  from  about  10^®  to 
10  -3 

4 X 10  cm  . Other  measurements  do  not  indicate  a minimum  although 

10-3 

densities  in  the  range  1.  5 to  5 x 10  cm  have  been  deduced  as  dis- 
cussed below.  Golomb  and  Good  suggest  the  minima  might  be  due  to  a 
very  complicated  wind  shear  and  turbulence  structure  in  the  90-130  km 
range  which  was  revealed  by  TMA  (tetramethylaluminum)  trails  released 
30  minutes  prior  to  their  experiment.  The  TMA  radiance  is  believed  to 
be  associated  with  the  atomic  oxygen  densities.  Thus  in  addition  to 
seasonal  and  latitudinal  variations,  local  motion  of  the  atmosphere  can  , 

apparently  contribute  significantly  to  variation  in  atomic  oxygen  density.  ^ 

i 

Atomic  oxygen  densities  have  been  deduced  from  satellite  meas-  q 

urements  of  the  5577  A airglow  [DG-73,  DG-74)  and  from  the  O-NO 
chemiluminescence  [Do-74].  The  photometer  data  clearly  demonstrate 
latitudinal  and  semiannual  variations.  But,  discounting  these  variations, 
the  deduced  oxygen  atom  densities  appear  to  be  significantly  smaller 
than  those  obtained  from  the  measurements  thus  far  considered.  The 
upper  limits  to  the  photometer  data  [DG-73,  Fig.  19]  are  plotted  on 
Fig.  2-1  and  the  range  of  all  measurements  is  shown  at  120  km;  these 
data  were  obtained  for  49  N.  There  is  some  uncertainty  in  the  value  of 
the  rate  constant  for  the  process  yielding  the  radiator,  0(^S),  but  to 
match  densities  observed  and  calculated  at  120  km  bv  others  (i.e.  , 

[O]  - 10  cm  ),  would  require  the  rate  constant  to  t)e  reduced  by  a 
factor  of  25.  (See  also  DG-74,  particularly  the  last  paragraph  of  the 
paper.) 

The  O-NO  chemiluminescence  measurements  ]Do-74|  are 
interpreted  as  leading  to  an  upper  limit  to  lN0]l0|.  Using  measured 


Y:  i 


values  of  [NO|,  the  O-atom  densities  have  been  computed;  a representa- 
tive profile  is  shown  on  Fip;.  2-1  [Do-74,  Fip.  1,  44  S,  lower  curve]. 
The  profile  is  ai?ain  consistently  lower  than  other  profiles  of  Fie;.  2-1. 

A very  recently  reported  measurement  |DB-74a|  in  the  altitude 
range  up  to  140  km  is  shown  on  Fig.  2-1.  These  were  nighttime  meas- 
urements at  57  N and  employed  a new  technique;  in  situ  measurement 
of  atomic  oxygen  by  absorption  and  resonance  fluorescence.  Absolute 
values  of  the  atomic  oxygen  density  were  determined  from  the  aljsorption 
experiments;  using  these  data,  the  resonance  fluorescence  measure- 
ments were  calibrated  and  were  extrapolated  to  oxygen  density  ranges 
where  the  absorption  technique  was  insensitive.  The  results  of  these 
measurements  [DB-74a]  are  in  reasonable  agreement  with  those  of 
Offermann  et  al.  |OD-73a,  OG-73  |. 

The  data  plotted  on  Fig.  2-1  are  for  a variety  of  times  in  the 
year.  Seasonal  variations  have  been  derived  from  data  analysis  by  sev- 
eral investigators,  at  least  [Sc-74a,  DG-73,  Of-74  |.  All  these  inves- 
tigations show  a minimum  in  the  atomic  oxygen  density  in  the  summer, 
between  June  and  September  and  also  demonstrate  yearly  variations  of 
about  a factor  of  three  in  atomic  oxygen  densities.  Of  the  measurements 
on  Fig.  2-1,  only  two  were  made  during  the  minimum:  Henderson's 
[He-74a]  in  August  and  Golomb  and  Good’s  |GG-72|  at  59  N also  in 
August.  The  remainder  were  made  in  March,  May,  September,  Octo- 
ber, and  November,  There  is  no  clear  seasonal  correlation  among  the 
measurements  shown  on  Fig.  2-1.  In  addition  to  the  summer  minimum, 
some  investigations  [Sc-74a,  DG-73  j show  a minimum  in  January; 

these  investigations  also  show,  inconsistentlv,  maximum  oxygen  atom 
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densities  of  2 or  3^10  cm  at  97  and  120  km,  respectiycly. 
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2.  2,  ? Model  Calculations 


Four  recent  model  calculations  are  considered  here  to  examine 
the  calculated  atomic  oxygen  profiles  for  altitudes  above  90  km.  No 
attempt  has  been  made  at  a thorough  comparison.  The  two  model  calcula- 
tions of  Shimazaki  and  Laird  [SL-70a,  SL-72]  and  Hunt  [Hu-73a]  give 
similar  results  as  indicated  on  Fig,  2-1.  Both  models  include  numer- 
ous chemical  reactions  (Hu-73a  has  the  greater  number,  150)  and  ac- 
count for  eddy  and  molecular  diffusion;  the  eddy  diffusion  coefficient 
profiles  (for  the  results  shown  on  Fig.  2-1)  are  also  similar.  Diffusive 
equilibrium  was  assumed  at  the  upper  boundary  (~150-160  km).  In 
spite  of  the  fact  that  these  calculations  were  canned  out  for  conditions 
at  the  equator  at  noon,  the  values  of  atomic  oxygen  density  are  in  rea- 
sonable agreement  with  the  measurements  of  Fig.  2-1  for  altitudes 
- 90  km  at  midlatitudes.  Shimazaki  and  Laird  [SL-72a]  have  also 
computed  O-atom  profiles  for  60  N,  summer-noon  conditions:  in  this 
calculation,  their  model  was  updated  for  cross-section  data,  tempera- 
ture profiles  and  the  eddy-diffusion  profile;  the  latter  profile  was  brought 
into  closer  agreement  with  obseiwations,  particularly  by  eliminating  the 
constant  value  of  the  eddy-diffusion  coefficient  above  100  km  and  replac- 
ing it  with  a coefficient  which  rapidly  declines  with  altitude  for  altitudes 
above  about  100-110  km.  The  oxygen  atom  profile  shown  in  Fig.  2-1 
now  reaches  a peak  just  above  90  km  and  at  higher  altitudes  declines 
more  rapidly  than  other  obseiwational  or  calculational  data  with  the  ex- 
ception of  the  data  of  Golomb  and  Good  |GG-72  |.  Although  calculations 
[SL-72a]  were  made  for  summer  and  winter  conditions,  onlv  the  summer 
data  are  plotted;  the  large  seasonal  difference  found  in  the  calculations 
[SL-72al  may  not  be  true  due  to  the  neglect  of  meridional  circulations. 
Note  that  differences  between  summer  and  winter  results  were  to  some 
extent  due  to  the  different  number  of  daylight  hours. 


The  model  calculations  by  George  et  al.  lGZ-72a|  are  shown 
on  Fig.  2-1  for  30  N,  noontime  with  a turbopause  at  94  km.  Their  major 
contribution  is  to  use  more  realistic  edclv  diffusion  coefficient  profiles 
and  to  examine  the  effect  of  changing  the  altitude  of  the  turbopause. 

Above  95  km,  agreement  is  not  good  with  the  observational  data  shown 
on  P'ig.  2-1  but  thev  would  probably  obtain  better  agreement  if  the  turbo- 
pause altitude  were  raised  as  I'an  be  deduced  from  their  results  for 
turbopauses  at  94  and  102  km  (with  the  exception  that  the  peak  U-density 
discrepancv  between  calculations  and  experiments  would  become  larger). 

In  the  model  calculations  of  Thomas  and  Bowman  lTB-72],  the 
assumption  of  diffusive  equilibrium  at  the  upper  boundarv  was  avoided 
for  the  species  O,  and  H2  by  allowing  fluxes  of  these  species  across 
the  upper  tioundary.  Their  eddv-diffusion  coefficient  profile  was  similar, 
qualitativelv.  to  that  of  Shimazaki  and  Laird  |S1  -72a  | for  summer  condi- 
tions but  had  a smaller  peak  value.  .\s  seen  on  Fig.  2-1.  their  calculated 
oxygen  atom  densities  are  clearly  lower  than  most  cither  plotted  data; 
the  conditions  of  the  calculations  were  45  N at  noon.  Thev  concluded 
that  a lower  eddy  diffusion  coefficient  would  vie  Id  better  .igreement  with 
the  preliminary  data  of  Henderson  |He-71  | but  their  eddv-diffusion  coeffi- 
cient was  already  smaller  than  that  of  George  et  al.  [GZ-72a]  whose  cal- 
culations are  in  better  agreement  with  the  final  data  of  Henderson  [lle- 
74a  [.  Clearly  a thorough  comparison  of  model  calculations  is  required 
to  aid  understanding  the  differences  in  model  calculations. 

A consistent  modelling  of  effects  of  diffusion  (molecular  and 
eddy)  on  neutral  composition  (O  included)  and  temperature  is  that  by 
Chandra  and  Sinha  |CS-74a|.  Thev  carried  out  a parameter  studv  and 
from  their  results  lCS-74a,  Figs.  2 and  4|  one  sees  that  the  m;iximum 
value  of  leddv  diffusion  coefficient)  and  the  altitude  of  the  turbo- 
pause, Z , have  little  effect  on  the  densities  below  100  km  (their  studv 
*^111 
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extends  to  85  km).  These  results  are  different  from  those  of  George 
et  al.  [GZ-72a)  where  a distinct  difference  in  the  atomic  oxygen  profile 
[GZ-72a,  Fig.  6]  occurs  for  various  Z^^^  at  altitudes  up  to  and  above 
100  km.  Furthermore  the  peak  O-densities  occur  at  about  90  km  in  the 
calculations  of  Chandra  and  Sinha  lCS-74al  for  Z = 105  to  120  km, 
but  at  105  or  87  km  in  the  calculations  of  George  et  al.  |GZ-72a|  for 
Z^^^  = 102  or  112  km,  respectively. 

Of  the  profiles  presented  by  Chandra  and  Sinha  [CS-74a|,  all 
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have  peak  densities  at  about  90  km  of  1.7  x 10  cm  . considerably 
greater  than  any  measurements  or  other  model  calculations  of  Fig.  2-1. 

2.  2.  3 Selection  of  Atomic  Oxygen  Profile 

In  selecting  an  atomic  oxygen  profile,  reliance  was  placed  on 
experimental- Results  for  altitudes  of  90  to  120  km.  At  each  selected 
altitude,  an  average  oxygen  density  was  computed  from  the  experimental 
values:  these  values  and  the  comjxited  averages  are  listed  in  Table  2-1. 
The  rocket  experiments  were  g-iven  a gi-eater  weight  than  the  satellite 
experiments  by  a factor  of  two  on  the  basis  of  the  more  direct  nature  of 
the  rocket  measurement  techniques.  Also,  in  accounting  for  tlie  satel- 
lite data,  not  only  were  the  upper  limit  oxygen  densities  derived  from 
photometer  data  [DG-73]  included  (see  above  discussion)  but  also  the 
lower  limit  oxygen  densities  of  these  photometer  data  [DG-73.  Fig.  19|. 

The  atomic  oxygen  profile  was  selected  liefore  the  pulilication 
of  the  extensive  measurements  of  Dickinson  et  al.  lDB-74al  and  these 
measurements  were  therefore  not  included  in  the  selection  process. 

Note  that  in  the  altitude  range  frf)ni  95  to  115  km,  the  data  of  Dickinson 
et  al.  jDB-74a|  deviate  by  no  more  than  15  percent  from  the  atomic 
oxygen  densities  of  the  selected  profile. 


'erage  Values  of  Oxygen  Atom  Densities  for  Altitudes  of  90  kni  and 
•eater  (derived  from  rocket  and  satellite  data). 


For  59°N  flight. 

For  22“N  flight. 

For  altitudes  > 100  km,  data  not  considered  as  a result  of  possible  local  atmospheric  motions. 


The  atomic  oxygen  pi'ofile  based  on  data  in  Table  2-1,  is  com- 


pared with  a model  of  Moe  lMo-73|  in  Table  2-2  and  Fig.  2-2  for  the 
altitude  range  90-120  km.  Moe  had  fewer  experimental  results  to  con- 
sider and  used  both  model  calculations  and  experimental  results  to 
arrive  at  his  model.  Both  model  profiles  have  peaks  at  about  97  km  but 
the  peak-density  of  the  distribution  of  Moe  is  about  twice  that  obtained 
in  the  present  review. 

2.  3 ATOMIC  OXYGEN  DENSITIES  BETWEEN  40  AND  90  KM 
2.3.1  Model  Calculations 

Below  90  km,  only  model  calculations  are  generally  available 
for  oxygen  atom  densities.  Those  calculations  which  take  into  account 
transport  phenomena,  extend  down  to  40  km  at  the  lowest  altitude  con- 
sidered. Below  40  km,  transport  processes  are  not  included  in  the 
available  calculations  which  yield  oxygen  atom  densities. 

Below  about  85  km,  a significant  diurnal  variation  in  atomic 
oxygen  densities  is  calculated.  For  the  present  models,  the  diurnal 
variation  will  be  simplified  by  selecting  one  model  for  noon  and  one  for 
midnight  conditions. 

Table  2-2.  Comparison  of  Atomic  Oxygen  Profile  of  Table  2-1 
with  Model  Profile  of  Moe  [Mo-73]. 


Atomic  Oxygen  Densities,  cm 
Altitude,  km  Table  2-1  I Moe 


2.  3.  1.  1 Daytime  Densities.  Between  40  and  90  km,  we  con- 
sidei'  the  five  recent  model  calculations  of  Hunt  [Hu-TSa],  Thomas  and 
Bowman  [TB-72],  Sliimazaki  and  Laird  [SL-70a,  SL-72,  SL-72a]  and 
George,  Zimmerman  and  Keneshea  lGZ-72a).  The  calculations  of 
George  et  al.  [GZ-72a]  show  that  profile  densities  are  relatively  insen- 
sitive between  80  and  90  km  to  the  tui-bopause  'Titude:  thus  in  this  range 
the  model  calculations  are  more  likely  to  agree.  This  is  shown  to  be 
the  case  in  Fig.  2-3  only  approximately.  At  80  and  85  km  all  calcula- 
tions apply  to  noon  and  there  is  no  correlation  with  altitude  of  the  spread 
in  calculated  densities.  Also,  no  simple  con-elation  with  gross  chemi- 
cal composition  seems  to  distinguish  the  calculations:  hydrogen-oxvgen 
systems  were  used  by  Thomas  and  Bowman  [TB-72]  and  George  et  al. 
[GZ-72a|,  a hydrogen-oxygen-nitrogen  system  by  Shimazaki  and  Laird 
[SL-70a,  SL-72,  SL-72a]  and  a hydrogen-oxygen-nitrogen-carbon  system 
by  Hunt  [Hu-73a|.  The  maximum  eddy  diffusion  coefficient  was  of  order 
10  cm  sec  in  the  calculations  of  George  et  al.  [GZ-72a|  and 
Shimazaki  and  Laird  [SL-72a  (and  SL-70a)]  whereas  in  the  calculations 

of  Thomas  and  Bowman  [TB-72]  and  Hunt  [Hu-73a],  the  m;tximum  was 
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of  order  10  cm  sec  . In  the  calculations  of  Shimazaki  and  Laird 
[SL-70a]  and  Hunt  [Hu-73a|,  the  eddy  diffusion  coefficients  were  given 
constant  values  above  100  km  and  this  is  not  reasonable  lSL-72a].  If 
we  choose  for  consideration  only  those  calculations  whicli  have  an  eddv 
diffusion  coefficient  profile  and  maximum  value  ci)nsistent  with  observa- 
tions [SL-72a],  then  the  calculations  of  Shimazaki  and  Laird  [SL-72a] 
and  George  et  al.  [GZ-72a]  are  selected.  That  this  is  a justifiable  pr' 
cedure  cannot  be  demonstrated  with  the  simple  comparisons  that  have 
been  presented;  a much  more  detailed  comparison  of  the  calculations  is 
required  and  perhaps  considerations  beyond  the  scope  of  anv  of  these 
calculations. 


Below  60  km,  the  calculations  of  George  et  al.  [GZ-72a]  show 
that  the  turbopause  altitude  has  no  effect  on  the  oxygen  atom  density. 

Thus  all  model  calculations  could  be  included  in  the  selection  process 
but  the  calculations  of  Thomas  and  Bowman  [TB-72|  and  Hunt  [Hu-73a] 
extend  only  to  60  km.  In  the  altitude  range  below  about  50  km,  the  cal- 
culations of  Nicolet  [Ni-71a]  are  important.  These  calculations  are 
based  on  photochemical  equilibrium  and  the  oxygen  atom  density  is  cal- 
culated by  using  an  ozone  profile  consistent  with  obseiwations  of  the 
ozone  densities. 

2.  3.  1.  2 Nighttime  Densities.  Of  the  recent  model  calcula- 
tions, only  those  of  Shimazaki  and  Laird  [SL-70a,  SL-72a|  and  Thomas 
and  Bowman  [TB-72]  yield  nighttime  profiles.  The  altitude  at  which 
the  calculated  oxygen  atom  densities  depart  from  the  daytime  values  is 
between  80  and  85  km.  In  all  cases,  the  oxygen  atom  density  is  negli- 
gibly small  by  60  km.  A recent  measurement  [DB-74a]  of  the  night- 
time atomic  oxygen  density  in  the  altitude  range  from  79.  6 to  90  km  is 
shown  on  Fig.  2-3.  The  departure  fi'om  the  daytime  values  between  80 
and  85  km  predicted  by  the  calculations  [SL-70a,  SL-72a,  TB-72]  is 
substantiated  by  the  measurements  [DB-74a|.  The  predicted  rapid  fall 

in  oxygen  density  is  also  confirmed  bv  the  measurements  which  show  the 
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density  to  be  less  than  3x10  cm  below  79.  6 km.  Additional  oxygen 
atom  densities  shown  on  Fig.  2-3  between  80  and  90  km  are  derived 
]RM-73|  from  OH  volume  emission  profiles. 

2.3.2  Selection  of  Atomic  Oxygen  Profile 

On  the  basis  of  the  calculations  of  Shimazaki  and  Laird  [SL-72a] 
for  summer  noon  at  60  N and  of  George  et  al.  [GZ-72a|  for  noon  at 
30  N,  an  average  daytime  oxygen  density  is  calculated  over  the  mutual 
range  of  altitudes  of  the  two  calculations,  i.  e.  , between  40  and  90  km. 
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In  the  ran^e  ol  altitudes  between  40  and  50  km,  the  photochemical  equi- 
librium calculations  of  Nicolet  [Ni-71a]  are  also  included.  Further- 
more, the  rocket  obseiwation  of  Golomb  and  Good  [GG-72|  at  85  km  for 
59  N is  included.  The  data  and  computed  averages  are  given  in 
Table  2-3. 

For  the  nighttime  oxygen  atom  profiles,  the  average  of  the 
values  of  the  calculations  of  Shimazaki  and  Laird  [SL-70a,  SL-72a| 
are  chosen  in  the  overlap  altitudes  of  these  calculations  as  given  in 
Table  2-3.  Below  75  km,  the  selected  oxygen  atom  densities  are  cal- 
culated on  the  basis  of  an  exponential  decline  from  the  atomic  oxvgen 
density  at  7 5 km  as  indicated  in  Section  2.  5 below.  The  calculations  of 
Thomas  and  Bowman  [TB-72]  are  for  pre-sunrise  and  post-sunset  con- 
ditions: however,  above  about  78  km,  the  profiles  for  pre-sunrise  and 
post-sunset  conditions  are  essentially  the  same  and  yield  lower  densities 
than  do  the  calculations  of  Shimazaki  and  Laird  |SL-70a,  SL-72a]  by 
up  to  three  orders  of  magnitude  at  78  km.  At  85  km,  these  calculations 
[TB-72,  SL-72a]  differ  by  a factor  of  2.  5. 

2.  4 ATOMIC  OXYGEN  DENSITIES  BELOW  40  KM 
2.4.1  Selection  of  Atomic  0.xyi:en  Profile 

Below  40  km,  the  calculated  oxygen  atom  density  continues  to 
decline  rapidlv.  We  relv  on  the  calculations  of  Nicolet  [Ni-71a[  to  an 
altitude  of  20  km:  see  Table  2-4a.  The  method  of  obtaining  these  den- 
sities has  been  mentioned  above.  The  calculations  are  for  an  overhead 
sun.  In  the  calculations,  no  significant  effect  of  hydrogenous  constitu- 
ents on  the  oxygen  atom  densities  was  found  below  50  km. 

Below  20  km,  the  calculations  of  Nicolet  lNi-72b|  were  ex- 
tended bv  assuming  equilibrium  between  photodissociation  of  and 


Table  2-4.  Calculated  Values  of  Atomic  Oxygen  Densities  Below  40  km. 


a.  20-40  km 


h.  km 

fO],  cm"^ 

Ref. 

35 

2.6(8) 

Ni-71a 

30 

4.5(7) 

Ni-71a 

25 

7.0(G) 

Ni-71a 

20 

1.0(C) 

Ni-71a 

b.  0-20  km 


h, 

km 

--J+ 

[O],  cm 

<[0]),  cm  ^ 

a 

b 

(a-fh)/2 

15 

2.61(5) 

1.2(5) 

1.9(5) 

10 

2.91(4) 

1.46(4) 

2.2(4) 

5 

7.7(3) 

4.74(3) 

6.  2(3) 

0 

6.6(0) 

4.9(0) 

5.  8(0) 

where 


[O]  = 


kLo. 


3X30"®  k'  J' 


1 + 


klO,]  [M] 


with 


k = rate  coefficient  for  reaction  O + ©2  ^ M -»  Oj  -*  M;  M = 02>  N2>  see  Ni-72b. 

k'  = rate  coefficient  for  reaction  NO  + Og  -*  NO2  ©2:  k = 10  exp(-1250/T),  see  Ni-72b. 

J = 0.  01  e.xp(-0.  38)  = 6.  84(-3)  for  sec  \ = 1,  see  Ni-72b.  N'02  + hi'  -•  NO  * O. 

j'  = 3.  7(-4);  Og  + ht'  (3000.4  < X < 11800A)  - O2  * 0(^P);  see  Le-72a. 


and 


h, 

km 

**■ 

tk 

0 

H 

1 

i 

[M],  cm  ^ 

I ^ 

K,  cm 

-r* 

sec 

— 1 

3 -r* 

cm  SGC 

[O3I. 

-3^ 

cm 

a 

b 

15 

217 

4.  05(18) 

1.  K-33) 

2.4(-33) 

3.  15(-15) 

2.60(12) 

10 
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recombination  of  O via  O + ©2  M.  The  ozone  densities  recommended 
by  Cadle  [Ca-74]  were  adopted  for  this  calculation  and  the  particle  den- 
sities were  from  the  U.  S.  Standard  Atmosphere  [US-621;  see 
Table  2-4b.  Levy  [Le-72a,  Le-73]  has  also  calculated  species  densi- 
ties in  the  troposphere  but  the  input  data  he  uses  are  not  the  most  reli- 
able data  available. 

2.  5 ATOMIC  OXYGEN  PROFILE  BETWEEN  0 AND  120  KM 

The  composite  oxygen  atom  densities  are  shown  in  Table  2-5 
and  Fig.  2-4  between  0 and  120  km  for  noon  and  midnight  conditions. 

Exti’apolation  of  the  profile  from  120  to  130  km  is  included  to  aid  fitting 
profiles  for  altitudes  above  120  km  to  the  profile  given  in  Table  2-5  and 
Fig.  2-4.  The  profile  is  for  midlatitudes  with  no  seasonal  or  latitudinal 
variations  taken  into  account. 
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Table  2-5.  Atomic  Oxypen  Profile  to  120  km  for  Noon  and 
Midnight  Conditions. 
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The  data  at  125  and  130  km  are  exl  rapol.itions  of  the  dat.i  at  lower 
altitudes  and  are  meant  to  b<’  used  to  • fittina  tlie  profile  lielow  120  km 
to  the  atomic  oxygen  profile  above  120  km. 

Noon  and  midniaht  values  are  identical  above  85  km. 

‘Based  on  [0|  [O  |„(.  exp|  - 1 . 33;75  - h)  | where  iO|„._  4.9(8)cm”^. 
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3.  DENSITY  PROFILE  OF  MOLECULAR  OXYGEN,  0„(a  A ) 


3.  1 


INTRODUCTION 


A midlatitude  profile  for  the  first,  excited  state  of  molecular 

oxygen,  0„(^A  ),  is  selected  for  altitudes  from  ground  level  to  beyond 
^ g 

140  km.  Fortunately,  measurements  have  been  made  between  30  and 

120  km  in  the  altitude  range  of  principal  interest  for  0„(^A  ).  At  other 

^ g 

altitudes,  as  well  as  the  edges  of  the  stated  range,  reliance  has  to  be 
placed  on  model  calculations.  The  measurements  of  02(^A^)  profiles 
have  been  mainly  made  at  latitudes  near  60  N so  that  latitudinal  or 
hemispherical  variations  in  the  profiles  have  not  been  established.  Also 
in  the  preliminary  stage  is  understanding  of  the  seasonal  variation  of  the 
0„(^A  ) densities. 

The  measured  profiles  of  02(^A^)  at  altitudes  above  approxi- 
mately 50  km  are  of  interest  also  in  connection  with  the  profiles  of 
ozone  at  these  higher  altitudes.  The  photodissociation  of  ozone  below 
90  km  is  the  principal  source  of  0„(  A ).  Since  measurements  of  ozone 

2 g 

densities  have  not  been  extensively  made  at  altitudes  in  excess  of  70  km, 
the  profile  of  0„(^A  ) between  70  and  90  km  provides  information  aimut 

2 g ^ 

the  ozone  density  in  this  altitude  range. 


3.2  MEASUREMENTS 

Measurements  of  the  (0,0)  band  at  1.27  um  from  the  radiative 
decay  of  02(^A^)  lEH-68|  were  made  between  30  and  110  km  on  a rocket 
flight  at  32  N latitude  in  October  1966.  These  measurements  are  shown 
on  Fig.  3-1;  a primary  peak  near  50  km  and  a secondary  peak  near 
85  km  are  obseiwed.  The  overhead  emission  on  this  flight  is  in  fair 
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agreement  with  an  earlier  balloon  measurement  [EL-67].  Above  80  km, 
the  detector  noise  for  the  rocket  flight  data  became  large  so  that  er- 
rors of  factors  of  2 and  5 at  90  and  110  km  were  possible.  The  rocket 
flight  measurements  of  Evans  et  al.  [EH-68|  were  later  normalized 
[EL-70a]  to  agree  with  the  balloon  obseiwations  [EL-69  |:  the  revised 
data  are  shown  on  Fig.  3-1  by  the  dashed  line  (note  that  the  shift  in  the 
secondary  peak  may  be  the  result  of  the  choice  of  data  represented  by 
the  solid  line  rather  than  by  the  dotted  line  of  Fig.  4b  of  EH-68).  Meas- 
urements at  the  higher  latitude  of  59  N by  Wood  et  al.  [WE-70],  pre- 
sumably using  the  same  procedures  as  Evans  et  al.  [EH-68],  have  been 
made  for  morning  and  afternoon  conditions  in  July  1967.  A shift  in  the 
primary  0„(^A^)  peak  density  to  higher  altitudes  (i.  e.  , 55  to  60  km)  is 
found  but  the  peak  densities  are  comparable  to  the  normalized  densities 
[EL-70a]  of  the  lower  latitude  (32  N)  data.  In  the  measurements  of 
Wood  et  al.  |WE-70],  evidence  for  a secondary  layer  between  80  and 
90  km  was  obtained.  Preliminary  data  of  Wood  [EL-70a]  for  evening 
measurements  are  shown  in  Fig.  3-1;  these  data  yield  lower  densities 


than  the  daytime  measurements  at  altitudes  below  about  90  km  but  the 

two-layer  structure  of  the  0„(/A  ) profile  is  again  found.  Companion 

^ g 

morning  measurements  [EL-70a]  between  70  and  100  km  are  not  sig- 


nificantlv  different  from  the  evening  measurements. 


Rocket  measurements  of  Megill  et  al.  [MH-70b|  between  70  and 
90  km  are  in  qualitative  agreement  with  the  normalized  values  of  02(^A^^) 
of  Evans  and  Llewellyn  [EL-70a]  liut  are  significantly  smaller  in  abso- 
lute values  between  76  and  84  km  while  having  satisfactory  aszreement 
otherwise,  as  shown  on  Fig.  3-1.  These  measurements  were  obtained 
at  59  N latitude  at  0100  UT. 
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The  secondary  0„(^A  ) layer  was  obseiwed  with  peak  intensity 
2 ^ • 

at  80  ±5  km  on  a rocket  flight  [HM-73^1  in  the  region  of  the  equator 

after  local  sunset;  the  peak  intensity  was  found  to  decay  with  time  after 

9-3  8 - 3 

sunset  and  yaried  from  about  1.  8 x 10  cm  to  8.  0 x 10  cm  during 
the  obseiwations.  An  apparent  latitudinal  variation  of  02^^A^)  emission 
was  detected  over  the  latitude  range  of  0 to  14  S. 

The  decay  of  both  layers  of  0„(^A  ) is  discussed  by  Vallance 

2 g 

Jones  [Va-73aj.  The  lower  layer  witn  peak  density  near  50  km  decays 
within  about  10  minutes  after  solar  illumination  ceases  whereas  the 
upper  layer  with  peak  density  near  85  km  decays  in  a characteristic  time 
of  about  60  minutes,  approximately  the  radiative  lifetime  of  02(^A^^). 

The  former,  rapid  decay  is  a result  of  quenching  at  the  lower  altitudes 
[EW-70J. 

Nighttime  measurements  [EL-72b,  BB-72c]  have  detected  only 
a single  0„(^A^)  layer  centei’ed  about  90  kno  Rocket  measurements  at 
59  N between  70  and  110  km  have  shown  a double  peaked  profile  as  ex- 
hibited on  Fig.  3-1;  these  investigators  [EL-72b|  consider  the  possi- 
bility that  a single,  broad  distribution  for  the  0„(^A  ) density  profile 
could  form  under  other  atmospheric  conditions.  By  contrast.  Bishop  et 
al.  lBB-72c]  found  only  a single  peaked  layer  between  70  and  110  km; 
these  investigators  considered  the  possible  enhancement  of  the  measured 
1.27  pm  radiation  by  auroral  phenomena  (although  a lack  of  enhancement 
is  indicated  by  aspects  of  the  observations).  The  enhancement  factor  is 
apparently  not  relevant  to  the  measurements  of  FJvans  et  al.  (f]L-72ijl. 

On  the  basis  of  ground-based  twilight  obseiwations  of  the  1.  27  ^'m 
radiation,  Evans  and  Llewellyn  [EL-72a;  see  also  EVV-70|  deduce  a 
seasonal  variation  of  the  02(^A^  ) layer  density  in  the  80-90  km  layer  for 
which  a minimum  occurs  in  mid-summer  and  a maximum  in  mid- winter. 
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The  variation  of  the  daytime  intensity  was  found  not  to  be  different  | EL- 
69].  A similar  seasonal  variation  has  been  reported  [VG-63|  from  ob- 
servations of  the  1.  58  iim  radiation  (the  (0, 1)  band  from  radiative  decay 
of  02(^A^)):  a maximum  in  intensity  was  found  around  January  and  a 
minimum  in  the  summer.  These  measurements  [VG-63]  were  evening 
twilight  measurements  at  52°N.  Note,  however,  that  a clear  seasonal 
variation  is  not  evident  in  the  compilation  of  twilight  and  day  airglow 
density  profiles  presented  by  Vallance  Jones  [Va-73a]  in  the  altitude 


range  of  80  to  90  km. 


For  the  lower  02(  A^)  layer,  no  seasonal  vari- 


ation is  evident  EL-72a  . 


Diurnal  measurements  of  02(  A^^)  densities  above  30  km  have 
been  made  by  Evans  et  al.  [EL-67,  EL-69J;  these  show  a relatively 
rapid  rise  and  decline  at  twilight  and  a slowly  varying  emission  during 
the  day  with  a peak  density  at  noon.  This  is  consistent  with  the  compari- 
sons of  Haslett  et  al.  [HM-69c]  in  which  the  integrated  intensity  meas- 
urements for  the  1.27  /im  radiation  of  several  groups  of  investigators 
for  restricted  inteiwals  of  solar  zenith  angle  are  combined. 


3.  3 


MODEL  CALCULATIONS 


Verv  simple  photo-equilibrium  calculations  by  Vlasov  [VT-69| 
and  bv  Evans  and  Llewellyn  lEL-70a]  are  in  good  agreement  with  meas- 
urements [EH-68,  EL-70a]  between  30  and  75  km.  These  calculations 
require  an  ozone  density  profile  as  input  but  otherwise  take  into  account 

only  the  photodissociation  of  ozone  and  the  de-excitation  of  0„(^A  ).  Note 

1 ^ It 

that  Evans  and  Llewellvn  [EL-70a]  assume  that  0„(  A ) is  formed  bv 

2 g 

photodissociation  of  ozone  at  wavelengths  greater  than  2660  A whereas 
Vlasov  [Vl-69]  may  have  accounted  for  photodissociation  at  wavelengths 
greater  than  2100  A.  Comparison  of  these  calculations  with  the  experi- 
mental data  is  shown  on  Fig.  3-2. 
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The  verv  extensive  calculations  of  Hunt  (llu-73al  vield  0„(  A ) 
profiles  with  density  values  sisinificantlv  smaller  than  obseiwations 
|EL-70a]  in  the  ran^e  75  to  85  km.  This  discrepancv  is  attributed 
jllu-73al  to  deficiencies  in  the  model  ozone  densities  which  are  appar- 
entlv  too  small  above  70  km.  Likewise,  the  calculations  of  Thomas  and 
Bowman  [TB-721  yield  values  of  density  sisinificantlv  smaller 

than  obsen-ed  values,  see  Fitt.  3-2,  and  this  discrepancy  can  apparently 
also  be  traced  to  small  values  of  the  calculated  ozone  density  above 
60  km.  However,  both  the  calculations  of  Hunt  |Hu-73aj  and  of  Thomas 
and  Bowman  [TB-72|  are  in  qualitative  agreement  with  obseiwations  and 
do  predict  the  occurrence  of  the  second  layer  with  peak  densitv 

at  85  km;  this  prediction  results  from  the  inclusion  of  hvdrogenous  spe- 
cies in  the  chemistrv  for  these  calculations.  The  calculations  of  Turco 
and  Sechrist  [TS-72b]  are  not  in  qualitative  agreement  with  the  obseiwa- 
tions  [EL-70a[.  Bv  contrast,  the  midnight  02(^A„.)  profile  of  Hunt 
jllu-73al  has  larger  density  values  in  the  altitude  range  between  85  and 
95  km  than  obseiwations  as  shown  on  Fig.  3-1. 

The  calculations  of  Shimazaki  and  Laird  [SL-72a]  for  60  N 
latitude  shown  on  Fig.  3-2  are  in  fair  agreement  with  the  measurements 
of  Wood  et  al.  1WF-70|  (see  Fig.  3-lb  The  peak  density  occurs  at 
60  km  in  the  calculation  and  at  55  and  60  km  in  the  afternoon  and  morn- 
ing measurements  of  WDod  et  al.  |V\E-70j;  since  the  measurements 
|WE-70l  were  in  Julv,  they  should  be  compared  with  the  summer  calcu- 
lations of  Shimazaki  and  Laird  [SL-72a]  (the  upper  points  in  the  densitv 
spread  as  shown  on  Fig.  3-2). 

Crutzen  et  al.  lC.T-71  | have  made  calculations  of  the  02'^A__) 

profile  from  0 to  70  km  using  recent  experimental  data.  Crutzen  ef  al. 

took  into  account  the  formation  of  at  wavelengths  liehnv  2660 

1 b 

as  well  as  formation  of  02!  A,^)  via  the  spin  forbidden,  ozone  photolvsis 
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at  wavelen^lhs  greater  than  3100  A;  the  latter  process  is  important  for 

production  at  low  altitudes  and  small  zenith  angles.  The  cal- 
culations  of  Crutzen  et  al.  [CJ-71]  were  made  for  zenith  angles  corre- 
sponding to  sec  6 values  of  4 and  1.4.  The  latter  corresponds  to  noon 
at  45  N in  spring  and  the  foi’mer  to  zenith  angles  for  the  rocket  experi- 
ments of  Evans  et  al.  [EH-68J.  The  calculations  [CJ-71)  and  experi- 
mental results  [EL-70a|  are  in  reasonable  agreement  above  40  km  as 
shown  on  Fig.  3-2. 

Calculations  of  the  diurnal  variation  of  densities  be- 

tween  40  and  80  km  have  been  made  bv  Vlasov  |Vl-70l;  these  are  in 
qualitative  agreement  with  obseiwations  [EL-67,  EL-69,  HM-69c]. 

3.4  SELECTION  OF  0„(^A^,)  PROFILE 

3.4.1  Daytime  Densities 

In  the  altitude  region  between  50  and  90  km,  the  corrected, 
observations  of  Evans  et  al.  [EH-68,  EL-70a]  are  selected  for  the  rec- 
commended  profile  at  noon  and  midlatitudes.  These  data  [EL-70a]  are 
midlatitude  data  (32  N)  but  were  obtained  in  the  afternoon.  However, 
calculations  by  Crutzen  et  al.  lCJ-71)  indicate  that  above  55  km,  the 
0„(^A  ) densities  do  not  differ  significantlv  over  the  range  of  zenith 
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angles  corresponding  to  the  time  between  noon  and  the  time  of  the  I'ocket 
flight  [EL-70a].  Also  note  that  the  rocket  flight  |EL-70a]  was  in  Octo- 
ber, roughly  midway  between  nuvximum  and  minimum  of  the  possible 
seasonal  e.xtremes  [EL-69,  EW-70,  EL-72a)  in 

Above  90  km,  the  Ogf^A  ) apparentlv  continues  to  decline:  the 
decline  has  not  been  measured  in  davtime  but  measurements  above  100  km 
have  been  made  at  night;  also  calculations  [Hu-73a|  have  been  extended 
to  160  km.  While  the  calculations  are  not  in  qualitative  agreement  with 
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observations  between  70  and  90  km,  at  90  km,  the  calculations  agree 
with  the  corrected  data  of  Evans  et  al,  |EH-68,  EL-70a|  and  are  in 
reasonable  agreement  with  nighttime  data  [EL-72b|;  the  calculations 
show  equality  of  day  and  night  0„(^A  ) densities  above  100  km  so  the 
previously  mentioned  agreement  is  further  evidence  that  the  calcula- 
tions above  90  km  may  be  adequate  for  the  purposes  of  profile  selection. 
Thus  above  100  km,  we  choose  the  data  from  calculations  by  Hunt 
[Hu-73a|.  At  100  km  the  nighttime  value  of  the  02(^A^)  density  [EL-72b| 
is  chosen  for  the  reference  cuiwe,  and  between  90  and  100  km  values  of 
the  densities  are  chosen  by  interpolation. 

Below  55  km,  the  calculations  of  Crutzen  et  al.  [CJ-71]  indi- 
cate that  the  noon  profiles  will  have  consistentlv  larger  values  of  ) 
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densities  than  those  obtained  in  the  afternoon  obseiwations  of  Evans  et 
al.  [EH-68|.  To  account  for  this,  the  calculations  of  Crutzen  et  al. 
ICJ-71 1 have  been  normalized  to  the  55  km  0„(^A  ) density  of  the  cor- 

2 g ■ 

rected  rocket  data  [EL-70a|  and  the  resulting  values  between  0 and 
50  km  chosen  as  reference  data;  these  data  are  plotted  on  Fig.  3-2 
where  they  may  be  compared  with  the  rocket  data. 

The  composite  data  for  the  noon,  midlatitude,  profile 

are  given  in  Table  3-1  and  shown  in  Fig.  3-3. 

3.4.2  Nighttime  Densities 

.■\bove  100  km  we  relv  on  the  calculations  of  Hunt  [Hu-73a| 
which  indicate  no  difference  in  0„(^A  ) densities  tietween  noon  and  mid- 
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night.  Below  100  km,  the  measurements  of  Evans  et  al.  |EL-72b|  are 
applicable  but  rather  than  carry  the  detailed  two-peak  profile,  we  will 
follow  the  calculated  profile  lHu-73a]  but  adjust  the  peak  to  near  the 
peak  values  observed  jEL-72b,  BB-72c|,  i.e.  . adjust  the  peak  nighttime 
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4.  DENSITY  PROFILE  OF  OZONE 


4.  1 INTRODUCTION 

The  densities  of  ozone  have  been  extensively  measui'ed  for  alti- 
tudes up  to  70  km.  These  measui'ements  demonstrate  important  vari- 
ations of  the  ozone  density  diurnally,  seasonally  and  with  latitude  and 
longitude  [Du-71,  Du-74,  PL-74al:  there  is  also  a significant  difference 
between  the  ozone  distribution  in  the  northern  and  southern  hemispheres 
[Du-70a].  At  lower  altitudes,  between  about  20  and  40  km,  vertical  and 
horizontal  transport  processes  dominate  the  distribution  while  at  greater 
altitudes,  photochemical  processes  are  more  important  [Ni-74l.  Re- 
views of  atmospheric  ozone  have  been  given  by  Diitsch  [Du-70a,  Du-71, 
Du-74  |.  While  fewer  measurements  have  been  made  above  70  km  than 
below,  a second  ozone  layer  has  been  obseiwed  above  80  km. 

From  the  available  data,  a midlatitude  ozone  profile  is  formu- 
lated for  noon  and  midnight  conditions  in  the  northern  hemisphere: 
seasonal  and  geographical  variations  are  not  included  in  tiie  model. 

4.  2 DAYTIME  DENSITIES 
4.  2.  1 Measurements 

From  the  daytime  obseiwations,  a midlatitude,  northern- 
hemisphere  model  of  the  ozone  profile  has  been  constructed  bv  Knieger 
and  Minzner  and  presented  prior  to  publication  t)y  Cadle  lCa-74|.  This 
model  applies  to  altitudes  up  to  72  km  for  a latitude  of  4 5 N.  In  the 
altitude  region  tjetween  0 and  26  km,  the  model  is  the  same  as  that  of 
the  U.  S.  Standard  Atmosphere  Supplements  jUS-66|.  Between  26  and 
34  km,  tjalloon  and  rocket  data  were  used  and  above  34  km,  rocket  data 
alone  were  used.  Data  from  onlv  17  of  37  rocket  flights  were  selected 
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for  use  in  the  model  on  the  basis  of  accuracy  of  measurements,  height 
resolution,  absence  of  abnormal  geophysical  conditions  and  deviation 
from  the  statistical  distribution  of  rocket  data.  Data  from  satellite  ob- 
servations based  on  backscattered  UV  radiation  were  used  to  latitudi- 
nally  adjust  the  mean  rocket  data.  The  data  employed  in  the  model  con- 
tain no  significant  longitudinal  gradients.  Because  of  the  time  at  which 
many  of  the  soundings  were  made,  the  data  for  altitudes  above  30  km 
are  biased  for  solar  conditions  in  the  years  1968  to  1970,  i.  e.  , for 
solar  minimum  conditions.  The  model  profile  is  plotted  in  Fig.  4-1. 


Between  72  and  80  km,  two  additional  measurements  lWS-68, 
RM-63)  provide  data  which  are  not  inconsistent  w’ith  the  model  atmos- 
phere profile  of  Krueger  and  Minzner;  these  data  are  shown  in  Fig.  4-1 
Also,  measurements  of  Krueger  as  given  in  Park  and  London  [PL-74a] 
from  altitudes  between  70  and  79  km  are  in  agreement  with  other  meas- 
urements shown  on  Fig.  4-1. 


Above  80  km,  there  are  tew,  direct  measurements  of  the  ozone 
densitv.  The  results  of  Rawcliffe  et  al.  lRM-631  extend  to  about  85  km 
and  indicate  a ledge  in  the  ozone  density  but  do  not  show  the  maximum 
in  ozone  densitv  which  recent  nighttime  [HR-72a,  HR-73,  SU-74]  or 
sunset  |MR-73l  measurements  reveal.  The  latter  measurements  of 
Miller  and  Ryder  |MR-73|  were  obtained  by  an  occultation  technique 
from  rocket  experiments  made  at  59  N latitude.  The  ascent  and  descent  ’ 

profiles  are  shown  on  Fig.  4-1;  the  ascent  profile  has  the  larger  peak 
value.  Those  profiles  have  peak  densities  between  86  and  89  km.  An- 
other ozone  profile  e.xtending  above  80  km  has  been  derived  |FL-70a]  ' 

from  a measured  0„(^A  ) profile  and  is  shown  on  Fig.  4-1  between  40  j 

and  90  km.  This  profile  generally  has  larger  densities  above  80  km  than 
either  of  the  profiles  obtained  from  other  e.xperiments  [RM-63,  MR-73] 
iiut  tielow  70  km,  is  in  good  agreement  with  the  model  profile  of  Kmeger  a 

and  Minznei'.  I 
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4.  2.  2 Model  Calculations 

Near  70  km,  the  calculations  of  Hunt  [Hu-73a|  are  in  close 
agreement  with  the  ozone  profile  of  Krueger  and  Minzner  [Ca-74|  and 
at  85  km  yield  a secondaiT  peak  in  ozone  density  which  is  in  qualitative 
agreement  with  the  ozone  profile  derived  from  measured  pro- 

files [EL-70a|.  The  comparison  of  these  profiles  can  be  made  from 
the  data  on  Fig.  4-1.  The  ozone  profile  derived  from  0„i^A  ) data 

2 g 

[EL-70a]  for  32  N has  larger  densities  than  the  calculations  of  Hunt 
[Hu-73aj;  the  latter,  however,  considers  his  densities  to  be  too  low 
above  70  km  based  on  a comparison  of  model  calculations  for  dawn  con- 
ditions with  satellite  measurements  lHR-72al.  Other  data  from  Evans 
and  Llewellyn  |EL-70a|,  for  a latitude  of  60  N,  yield  profiles  in  the  75 
to  90  km  range  which  exhibit  a considerable  variation  of  densities:  but 
the  lowest  density  at  80  km  is  still  twice  as  large  as  the  computed  den- 
sity of  Hunt  [Hu-73a]. 

The  calculations  of  Thomas  and  Bowman  [TB-72]  extend  from 
60  km  to  125  km;  between  60  and  75  km,  these  calculations  are  in  good 
agreement  with  the  model  ozone  profile  lCa-74|.  Above  75  km.  the 
calculations  agree  moderately  well  with  the  calculations  of  Hunt  |Hu-73a]. 
being  somewhat  lower  between  7 5 and  100  km  and  slightly  higher  between 
100  and  125  km. 

The  calculations  bv  Shimazaki  and  Laird  [SL-72a]  for  60  N for 
summer  and  winter  conditions  yield  densities  greater  than  the  ozone 
model  profile  for  altitudes  70  km  but  practicallv  span  the  other  cal- 
culations [Hu-73a,  TB-72|  and  obsen-ations  |WS-68,  RM-63.  EL-70a, 
MR-73]  at  altitudes  between  70  and  91  km.  The  previous  calculations 
of  Shimazaki  and  Laird  [SL-70a|  forO  latitude  yield  S'i gnificantlv 
smaller  densities  than  those  shown  on  Fig.  4-1. 
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Park  and  London  [PL-74a|  have  also  calculated  ozone  profiles 
up  to  100  km;  the  profile  in  the  altitude  range  between  7 5 and  100  km 
does  not  qualitatively  agree  with  the  measured  [EL-70a,  MR-73  | oi 
calculated  [Hu-73a,  TB-72]  profiles  as  shown  on  Fig.  4-1. 

4.2.3  Selection  of  Ozone  Profile 

For  altitudes  between  0 and  75  km,  the  model  ozone  profile 
from  Cadle  [Ca-74]  is  selected.  Between  70  and  80  km,  several  of  the 
rocket  and  satellite  expei'iments  [RM-63,  WS-68,  Krueger,  1968  — see 
PL-74al  are  in  rough  agreement  and  also  roughly  agree  with  an  extrap- 
olation of  the  model  ozone  profile  to  about  80  km.  The  ozone  profile 
derived  from  measurements  of  the  02(^^g.)  emission  [FL-70a]  has 
larger  0„  densities  than  the  other  measurements  between  75  and  80  km, 

O 

yet  does  exhibit  the  maximum  at  about  85  km  which  the  comprehensive 
calculations  [Hu-73a,  TB-72  1 predict.  (As  discussed  below  - Sec- 
tion 4.3.  1 — the  nighttime  measurements  [HR-73,  SU-74|  also  reveal 
a maximum  in  ozone  density  near  85  km.)  However,  these  calculations 
[Hu-73a,  TB-72J  predict  lower  densities  at  80  km  than  the  set  of  ob- 
servRions  [RM-63,  WS-68,  Krueger,  1968-  see  PL-74alat  80  km.  As  a 
compromise,  the  average  of  the  rocket  and  satellite  obseiwations  be- 
tween 75  and  80  km  will  be  used  and  between  80  and  90  km  the  profile 
,vill  be  the  average  of  the  indirect  obseiwations  lEL-70al  and  the  calcu- 
lations [Hu-73a,  TB-72,  SL-72al.  Above  90  km,  the  average  of  the 
calculations  will  be  chosen  [llu-73a,  TB-72  |.  The  input  data  and  calcu- 
lated averages  are  shown  on  Table  4-1,  (The  measurements  of  Miller 
and  Ryder  (MR-731  were  not  included  in  the  selection  of  the  ozone  pro- 
file as  they  were  unknown  to  us  at  the  time;  bv  including  the  average  of 
the  ascent  and  descent  mea.su rements,  the  average  ozone  densitv  would 
change  l>v  less  than  15  l)etween  75  and  90  km  and  bv  -33  at  90  km.) 
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Table  4-1.  Calculation  ol  Average  Ozone  Densities 
Between  7 5 and  125  km. 


IO3  |.  cm 

3 

Kr;S 

RM-63 

WS-68  PL-74a  EL-70a 

Hu-73a  TB-72 

SL-72a'^ 

1.3(8)* 

1.7(8)  1.0(8) 

1.0(8)'’’ 

4.4(7)"  1.5(8) 

2.65(7)  1.95(7) 

2.  2(8) 

5.6(7)  3.3(7) 

1.8(8) 

4.  0(7) 

1.6(7)  1.0(7) 

5.  2(7) 

2.9(6)  1.85(6) 

2.5(5)  2.4(5) 

2.  5(4)  i 4. 2(4) 

i 

1 

2.8(3)  4.6(3) 

2.9(2)  4.6(2) 

3.7(1)  4.6(1) 

1.0(1)  4.6(0) 

Work  of  Krueger  reported  in  reference  PL-74a. 
Average  value. 

Interpolated  value. 

E.xtrapolated  value. 

Average  of  summer  and  winter  values. 


^[03!) 

-3 

cm 

1.3(8) 
6.  8(7) 


The  selected  data  are  plotted  on  Fig.  4-2  and  listed  in  Table 
4-2.  The  secondary  layer  in  the  density  centered  about  85  km  has 
a strong,  seasonal  variation  [EL-72a|,  as  deduced  from  twilight  meas- 
urements of  the  0„(^A  ) emission;  this  layer  has  a maximum  ozone 
density  in  midwinter  and  a minimum  density  in  midsummer.  The  mean 

annual  concentration  in  a 10-km  layer,  as  deduced  from  the  obseiwed 
1 7-3 

0„(  A ) emission,  is  8 x 10  cm  . whereas  the  mean  concentration  of 

the  selected  profile  shown  in  Fig.  4-2  in  a 10-km  layer  centered  at 
7-3 

85  km  is  7 X 10  cm  . Thus  the  selected  profile  in  this  10-km  layer 
represents  a type  of  mean  annual  profile.  Note,  however,  that  Hayes 
and  Ruble  [HR-73]  found,  in  the  equatorial  regions,  no  apparent  sea- 
sonal pattern  in  the  ozone  profiles  between  50  and  100  km  for  nighttime 
measurements. 

4.  2.  4 Late  Reports  on  Ozone 

Recently  [Tr-75]  a daytime  measurement  of  the  ozone  density 
in  the  altitude  region  between  90  and  110  km  has  been  successfully  made 
for  the  first  time  by  use  of  a cryogenic,  neutral  gas,  mass  spectrometer. 
Preliminary  results  are  shown  on  Fig.  4-1;  these  data  are  upper  limit 
values  and  could  be  smaller  by  60  o at  90  km  but  only  by  5 ■ at  110  km. 
Above  100  km,  the  measured  densities  are  significantly  larger  than  the 
calculated  densities  of  Hunt  [Hu-73a]  and  Thomas  and  Bowman  [TB-721. 
These  ozone  measurements  ] Tr-75  | are  put  on  an  absolute  basis  by- 
using  the  CIILA.-72  density  of  argon  at  90  km  as  a normalizing  point. 

4.3  NIGHTTIME  DENSITIES 
4.  3.  1 Measurements 

The  nighttime  ozone  density  only  varies  significantly  from  the 
daytime  ozone  density  at  altitudes  above  about  55  km,  according  to 
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Table  4-2.  Profile  of  Ozone  for  Noon  and  Midnight  Conditions. 


theoretical  expectations  [Ca-74  | and  accordino;  to  observations  |CH-66a, 
Re-68c,  Hi-71a,  CH-72j.  This  is  demonstrated  on  Fig.  4-3  where  the 
ozone  profile  selected  as  the  daytime  reference  is  drawn  as  a smooth 
cuiwe.  At  altitudes  below  about  55  km,  there  is  excellent  agreement 
between  the  nighttime  measurements  (CH-66a,  Re-68c,  Hi-71a, 

CH-72|  and  the  daytime  ozone  profile  with  the  exceiJtion  of  satellite 
measurements  [HR-73  j;  above  55  km,  the  nighttime  obseiwations  yield 
larger  densities.  In  the  vicinity  of  65  km,  there  is  a large  spread  in 
the  experimental  values.  The  rocket  firings  [CH-66a,  Re-68c,  Hi-71a, 
CH-72),  in  the  latitude  range  30-38°N,  yielded  data  up  to  about  70  km. 
These  experiments  used  as  a radiation  source  lunar  UV  radiation 
[CH-66a,  CH-72]  and  airglow  [Re-68c|  or  employed  a chemilumines- 
cent ozone  sonde  [Hi-71aJ.  The  satellite  obseiwations  [HR-73  |,  between 
50  and  80  km,  have  ozone  densities  smaller  than  the  daytime  profile  and 
this  is  in  contrast  to  all  other  measurements  as  well  as  calculations 
(see  below).  However,  these  measurements  e.xtend  to  100  km  and  do 
exhibit  the  secondaiw  ozone  peak  at  85  km.  Recent  nighttime  rocket  ex- 
periments [SU-74  I in  which  emission  measurements  of  ozone  at  9.6  pm 
have  been  made,  indicate  a peak  ozone  density  at  about  86  km. 

4.  3.  2 Model  Calculations 

The  calculations  of  Hunt  |llu-73a|,  Thomas  and  Bowman  (TB- 
72],  Shimazaki  and  Laird  [SL-70a,  SL-72a]  and  Park  and  London 
|PL-74a|  are  shown  on  Fig.  4-3.  Between  55  and  75  km.  the  calculated 
nighttime  ozone  densities  exceed  the  davtime  values  of  the  model  ozone 
profile;  between  75  and  85  km,  the  comjxited  densities  lHu-73a.  TB-72] 
are  smaller  than  the  davtime  profile  densities  and  at  altitudes  aiiove  85  km, 
all  nighttime  calculated  ozone  densities  are  larger  than  the  davtime  pro- 
file densities.  The  calculations  span  the  latitude  range  from  zero  to  60  . 
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The  total  ozone  content  of  the  atmosphere  increases  with  increasini^  lati- 
tude [Du-74  J,  however  this  is  not  generally  tnie  for  the  ozone  content 
at  specific  altitudes.  Between  about  30  and  50  km,  the  ozone  density  de- 
creases with  increasing  latitude  and  betw'een  60  and  80  km  this  depend- 
ence is  found  for  January  but  the  reverse  dependence  is  found  for  Julv 
[PL-74aJ.  Between  50  and  60  km  there  is  little  dependence  on  the  lati- 
tude. See  Fig.  4-4.  The  calculations  of  Shimazaki  and  Laird  [SL-70a, 
SL-72a)  are  in  qualitative  agreement  with  these  trends  for  altitudes  be- 
tween 40  and  70  km;  at  higher  altitudes  an  increase  in  the  ozone  density 
with  latitude  is  predicted  for  summer  and  winter  seasons.  The  calcula- 
tions of  Hunt  [Hu-73a|  at  0 in  comparison  with  those  of  Shimazaki  and 
Laird  [SL-72a]  at  60  , results  in  a prediction  of  an  increase  in  ozone 
density  with  latitude. 

The  diurnal  variation  of  ozone  has  been  calculated  [Hu- 73a, 
TB-72,  PL-74a);  Park  and  London  [PL-74a]  point  out  that  diurnal  asym- 
metry results  from  inclusion  of  hydrogenous  compounds  in  the  chemical 
system. 

4.3.3  Selection  of  Ozone  Profile 

For  altitudes  of  55  km  and  smaller,  the  nighttime  ozone  profile 
wall  be  taken  to  be  identical  to  the  da\dime  profile.  This  is  supported 
by  the  e.xcellent  agreement  between  the  selected  daytime  ozone  profile 
and  the  nighttime  observations.  Between  55  and  75  km,  reliance  will  be 
placed  on  the  obseiwations  and  the  average  of  the  obseiwed  values  will 
be  used  as  the  selected  ozone  nighttime  profile  (with  the  exception  of 
satellite  obseiwations  [HR-73]  which  are  significantly  different  from  the 
total  sets  of  calculations  and  remaining  observations!.  .>\t  and  above 
75  km,  the  average  of  the  recent  calculations  [TB-72,  SL-72a,  Hu-73a, 
PL-74aJ  will  be  used.  In  the  case  of  the  calculations  of  Shimazaki  and 
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Laird  [SL-72a]  the  average  of  the  summer  and  winter  calculations  is 
used;  this  has  an  effect  of  converting  the  calculations  [SL-72a|  to  a 
lower,  effective  latitude  than  the  60  N latitude  of  the  calculations  since 
the  ozone  density  in  summer  (upper  point  on  Fig.  4-3  for  SL-72a)  will 
be  lower  at  midlatitudes  than  at  60  N [PL-74aJ  and  the  ozone  density  in 
winter  (lower  point  on  Fig.  4-3  for  SL-72a)  will  be  higher  at  midlati- 
tudes than  at  60  N [PL-74a],  at  least  between  60  and  80  km.  Further- 
more, at  midlatitudes  (say  4 5 N)  the  seasonal  variation  in  total  ozone 
density  will  peak  near  summer  and  show  a minimum  near  winter  [Du-74], 
so  that  an  average  is  appropriate  for  a mean,  midlatitude  ozone  model. 

The  average  values  are  calculated  from  data  in  Table  4-3  for 
h 2:  60  km.  The  selected  nighttime  profile  is  shown  in  Fig.  4-2. 
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correcting  tins  error. 


5.  DENSITY  PROFILE  OF  NITRIC  OXIDE 


5.  1 INTRODUCTION 

A midlatitude  profile  for  nitric  oxide  is  selected  in  the  altitude 
range  from  ground  level  to  240  km.  While  measurements  of  nitric 
oxide  have  been  made  in  a large  fraction  of  this  altitude  range,  signifi- 
cant discrepancies  and  unceilainties  remain  among  the  data.  Simi- 
larly, agreement  among  calculated  nitric  oxide  profiles  is  lacking. 
However,  in  at  least  one  instance  of  apparent  discrepance  between  ex- 
perimental results,  the  difference  in  season  and  solar  activity  at  the 
time  of  the  measurements  has  been  shown,  by  model  calculations  [OJ-74], 
to  be  important  in  accounting  foj'  the  discrepancy.  Discussions  of  the 
effect  of  solar  activity  and  season  on  nitric  oxide  densitv  are  presented 
in  Sections  5.  2.  3 and  5.  4.  2.  Additional  measurements  and  model  calcu- 
lations are  necessary  for  establishing  reliable  nitric  oxide  atmospheric 
densities. 

5,  2 DAYTIME  NITRIC  OXIDE  DENSITIES  ABOVE  120  KM 
5.  2.  1 Measurements 

There  have  been  two  measurements  at  altitudes  above  120  km 
from  which  daytime  NO  densities  are  obtained;  a rocket  experiment  at 
32  N |FT-74]  in  which  the  twilight  emission  from  the  >-  and  6-liands 
of  nitric  oxide  were  obseiwed  and  a satellite  experiment  |Kr-70a]  in 
which  the  y-  and  6-bands  were  also  measured.  Near  215  km  where 
NO  densities  have  been  derived  from  the  two  measurements,  there  is 
about  an  order  of  magnitude  scatter  in  the  densities  as  shown  on  Fig.  5-1. 
This  scatter  is  discussed  in  Section  5.  2,3  below:  the  nighttime  measure- 
ments shown  on  Fig,  5-1  arc  discussed  in  Section  5.5  below. 
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5.  2.  2 Model  Calculations 


A number  of  model  calculations  have  been  made  [St-71,  SH-70, 
St-72a,  Hu-73a,  SL-70a,  SL-72a,  Ru-73a,  Ni-70,  O.T-74]  for  altitudes 
above  120  km.  In  the  calculations  of  Strobel  [St-72a|  shown  on  Fi".  5-1, 
boundary  conditions  were  carefullv  treated  but  a constant  eddv  diffusion 
coefticient  [St-71 1 was  used;  the  latter  is  not  in  accord  with  obseiwa- 
tional  data.  The  calculation  of  Hunt  [Hu-73a]  employed  an  eddv  diffusion 
coefficient  which  varied  with  altitude;  here  the  boundaiy  conditions  at 
the  upper  boundary  were  different  from  those  of  Strobel.  Hunt  [Hu-73a] 
used  a diffusive  equilibrium  condition  where  Strobel  [St-72a]  used  a 
flu.x  condition. 

The  calculations  of  Rusch  [Ru-73a]  yielded  larp;er  densities 
than  obseiwed  or  calculated  bv  others  and  the  calculations  of  Shimazaki 
and  Laird  [SL-70a,  SL-72aj  are  siy;nificantlv  smaller  since  ionic  reac- 
tions were  not  included.  The  calculations  of  Nicolet  lNi-70]  do  not 
include  transpoi-t  effects.  These  calculations  [Ni-70,  Ru-73a.  SL-70a. 
SL-72a]  were  omitted  from  the  comparisons  of  Fir;.  5-1. 

Recentlv,  calculations  [OJ-74|  have  been  made,  based  on  a 
model  similar  to  previous  models  [SM-70,  St-71,  St-71d]  but  with  an 
improved  description  of  chemical  processes  for  higher  altitudes  and  an 
altitude-dependent,  eddy  diffusion  coefficient.  These  calculations  were 
specifically  made  for  the  conditions  on  the  days  of  the  rocket  flights  of 
Meira  |Me-71[  and  Tisone  [Ti-73|;  the  noon  profiles  for  altitudes  aljove 
100  km  are  shown  on  Fig.  5-1.  The  noon  profiles  of  Oran  et  al.  |OJ-74| 
are  in  agreement  with  the  calculations  of  Hunt  |Hu-73a|  for  altitudes  up 
to  160  km.  At  higher  altitudes,  a significant  difference  in  nitric  o.xide 
densitv  for  the  conditions  of  the  two  flights  [Me-71,  Ti-73|  is  predicted 
bv  Oran  et  al.  |OJ-74[;  this  difference  is  further  considered  in 
Section  5.  2.  3. 
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5.  2.  3 Further  Analysis  of  Measurements  and  Calculations 

A large  disparity  among  measurements  at  215  km  has  been 
noted.  Also,  mention  has  been  made  of  recent  calculations  |OJ-74  ) in 
which  the  influences  of  solar  activity  and  seasonal  variations  on  the 
nitric  oxide  density  are  considered.  These  influences  will  be  further 
examined  here  before  attempting  to  select  a nitric  oxide  profile. 

In  the  vicinity  of  the  altitude  of  215  km,  there  are  two  da\"time 
measurements  from  which  nitric  oxide  densities  have  been  derived 
[Kr-70a,  FT-74j.  In  addition,  recent  calculations  [OJ-74|,  made  spe- 
cifically for  the  conditions  of  two  other  measurements  at  lower  altitudes 
[Me-71,  Ti-73  ],  have  been  extended  to  21 5 km.  By  comparing  the 
measured  and  calculated  nitric  oxide  densities  with  the  10.7  cm  solar 
flux  for  the  days  or  measurement  period  of  the  four  related  experiments 
[FT-74,  Kr-70a,  Me-71,  Ti-73],  a correlation  is  found  as  shown  bv  the 
data  of  Table  5-1  and  their  graphical  representation:  note  circles  on 
Fig.  5-2.  The  cuiwe  on  Fig.  5-2  will  be  discussed  later.  The  datum 
of  Feldman  and  Takacs  [FT-74)  is  reported  as  an  upper  limit  and  is 
shown  as  such  on  Fig.  5-2.  Note  that  the  calculations  of  Oran  et  al. 
[OJ-74  j are  included  in  Fig.  5-2  on  the  basis  that  they  are  assumed  to 
predict  reasonably  well  the  nitric  oxide  densities  at  high  altitudes  and 
not  on  the  basis  that  the  solar  intensities  used  in  the  calculation  10J-74] 
are  corrected,  in  some  manner,  to  correspond  to  values  of  the  solar 
flux  associated  with  the  calculated  nitric  oxide  densities  shown  on 
Fig.  5-2.  The  text  [OJ-74  | is  not  explicit  on  this  point  but  apparentlv 
no  such  correction  was  made.  The  rationale  for  the  relation  between 
nitric  oxide  density  and  the  solar  flux  is  as  follows.  According  to 
Jaccliia  [.Ia-71c|,  the  exospheric  temperature,  and  consequentlv  the 
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Table  5-1.  Nitric  Oxide  Densities  at  215  km  and 
the  Decimetric  Solar  Flux 


Solar  10.  7 cm  Flux, 

1 

) 

Day  or 
Period 

^10.  7 

_99  _2  -1 

10  W m Hz  ^ 

Ref. 

[NO) 

-3 

cm 

Ref. 

Jan  4,  1974 

75 

a 

^ 7.  0(5) 

FT-74 

June  1968 

142 

b 

5.  0(6) 

Kr-70a 

Jan  31,  1969 

130 

1 2.  9(6) 

1 

Me-71 

Feb  6,  1969 

142.  5 

d 

Me-71 

May  26,  1971 

91 

e 

1.  5(6) 

Ti-73  i 

^a.  "Solar-Geophysical  Data,"  354,  Part  I,  7,  Feb.  1974. 


b.  Ibid,  lER-FB-287,  7,  July  1968. 

c.  Ibid,  IER-FB-294,  7,  Feb.  1969. 

d.  Ibid,  IER-FB-295,  7,  Mar.  1969. 

e.  Ii)id,  322,  Part  I,  7,  June  1971,  U.  S.  Department  of  Commerce 
(Boulder,  Colorado  80302). 


molecular  oxygen  density  at  high  altitudes,  is  proportional  to  the  10.7 
cm  solar  flux.  But  the  nitric  o.xide  density  is  also  proportional  to  the 
molecular  oxygen  density  lOJ-74|;  therefore  the  nitric  oxide  density  is 
proportional  to  the  decimetric  solar  flux.  The  basis  of  the  dependence 
of  the  nitric  oxide  density  on  the  molecular  oxvgen  density  is  the  follow- 
ing reaction  couple; 


4 

N(^S) 

^^2 

> NO  + 0 

ia) 

4 

N(  S) 

• NO 

. N'2  • O 

lb) 

An  increase  in  moleculai'  oxvgen  increases  the  formation  rate  of  NO  via 
reaction  1 and  reduces  the  loss  rate  of  NO  in  reaction  b as  a result  of 
the  increased  loss  rate  of  N(^F)  in  reaction  a. 
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10  Altitude  ^ 215  km 


s / 

y 

/ 


i 1 i- 


-22  , -2  -1 

Decimclric  Solar  Flux,  10  W m Hz 


Fic:.  5-2.  Nitric  Oxide  Densities  vs  Decimctric  Solar  Flux. 


At  lower  altitudes,  between  120  and  160  km,  the  available 
measurements  lOG-73,  Kr-70al  and  calculations  [Hu-73a,  OJ-741  are 
in  much  better  agreement  than  at  215  km  so  that  the  strong  dependence 
of  the  nitric  o.xide  density  on  the  10.7  cm  solar  flux  is  not  e.xpected  at 
these  lower  altitudes.  This  is  consistent  with  the  fact  that  the  variation 
in  the  molecular  oxygen  densitv  with  variation  in  the  exospheric  tem- 
peratui’e  (and  thus  solar  flux)  is  smaller,  bv  factors  of  220  to  15,  in  the 
altitude  range  between  120  and  160  km  than  at  220  km,  according  to  the 
Jacchia  atmosphere  [Ja-71c|.  Consequentlv,  at  120  km,  the  assumption 
will  be  made  that  the  nitric  oxide  densitv  is  independent  of  the  decimetric 
solar  flux  (see,  however.  Section  5.4.2  below). 


64 


r 


Another  possible  influence  on  the  oxygen  density,  and  hence  on 
the  nitric  oxide  density,  is  seasonal.  Roble  and  Norton  [RN-72al  find 
a seasonal  variation  in  the  oxygen  densities  for  altitudes  in  the  range 
100  to  120  km.  The  measurements  are  from  solar  occultations  per- 
formed over  a nine  month  period  and  show  a minimum  molecular  oxygen 
density  in  winter  and  a maximum  density  in  the  summer  about  two  to 
three  times  larger  than  the  winter  minimum.  Previous  collations  of 
data  |Ch-70b)  have  shown  no  seasonal  variation  at  midlatitudes.  At 
altitudes  above  120  km,  apparently  no  seasonal  variation  has  been 
established. 


5.  2.  4 Selection  of  Nitric  Oxide  Profile 


Since  the  variation  of  nitric  oxide  density  with  solar  activity 
appears  to  be  strong  and  to  be  useful  in  accounting  for  the  disparity  in 
measured  and  calculated  densities,  the  selected  profile  will  include  the 
dependence  of  density  on  solar  activity.  .*\.ny  possible  seasonal  variation 
of  nitric  oxide  density  will  l)e  neglected.  At  215  km,  the  relationship 
between  density  and  solar  activity  will  be  roughly  appro.ximated  by  use 
of  the  following  equation; 


ln[NO]2i5 


6,  08  F^q  r. 

10.  22  + E 

5X10  - 


(1) 


-22  -2  -1 

where  F^q  ^ is  the  solar  activity  in  units  of  10  W m Hz  . The 

dashed  cutwe  on  Fig.  5-2  is  i)ased  on  Eq.  (1).  At  120  km,  the  nitiic 

T 3 

oxide  density  will  be  fixed  at  3 > 10  cm"  based  on  obseiwation  lKr-70al. 
To  establish  the  nitric  oxide  profile  at  and  aiiove  120  km,  the  slopes 
(-  -d  InjNOl  dh)  will  be  determined  bv  the  relations: 
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s = s 


ln;iNO]2j5/lNOl^.| 

215  - h 

r 


h = 120  km  (2) 

120  km  < h ^ (3) 

h ^ h (4) 


where  S is  the  slope  at  120  km,  the  variable  slope  between  120  and 

h km,  and  S the  constant  slope  at  and  above  h km.  Bv  using  the 

definition  of  slope  and  Eq.  (3),  the  nitric  oxide  densitv  between  120  and 

h km  is  found  to  be 
r 

2 

InlN-Ol  = ln[NO|,2g-  (h-  120)S  - | '|h  ‘ (=> 

and  at  h^  km,  the  density  is  thus 

(h  - 120) 

InLNOl^,  = ln[NO]^20  ■ * ^c^ 

At  altitudes  above  h 

r 

In'iNOl  = InlNOl  - S (h  - h ) iT) 

L j ' 'r  c r 

By  using  Eqs.  (4)  and  (6),  the  slope,  S^.,  becomes 

(2  1n.;iNOl2,5  lNOl,,„l-lh,.-  120)S  I 
" ■ I 310  - j 

Profiles  for  nitric  oxide  at  altitudes  above  120  km  can  be  calculated  bv 
using  Eqs.  (5)-(8)  togetlier  with  the  selected  nitric  oxide  densitv  at 
120  km,  the  nitric  oxide  densitv  at  215  km  as  calculated  with  Eq.  '!). 
and  the  slope  S at  120  km  as  determined  from  the  nitric  oxide  profile 
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fur  120  km  and  lower  altitudes.  (See  Section  5,  3.  3.  ) The  profiles  of 
nitric  oxide  above  120  km  are  illustrated  in  Fig.  5-3  where  h^,  = 125  km. 

5.3  DAYTIME  NITRIC  OXIDE  DENSITIES  BELOW  120  KM 
5.  3.  1 Measurements 

A number  of  obseiwations  have  been  made  [Ba-66d,  Pe-69a. 
Me-71,  Ru-73a,  Ti-73]  at  altitudes  above  60  km  by  measurement  of  the 

radiation  from  the  NO  > -bands.  For  all  but  one  investigation  |Pe-69a]. 

7 8-3 

the  NO  densities  lie  between  10  and  10  cm  ; the  results  of  Pearce 
[Pe-69a]  as  can  be  seen  on  Fig.  5-4  are  more  than  an  order  of  magni- 
tude larger  than  the  other  observations.  The  profiles  deduced  by  Pearce 
[Pe-69aj  have  been  criticized  by  Meira  |Me-7l]  on  the  basis  of  inade- 
quate correction  for  Rayleigh  scattering  and  ozone  absorption  contributions. 

From  the  measurements  of  Rusch  [Ru-73a]  onlv  upper  and 
lower  limits  are  ol)tained;  these  are  shown  on  Fig.  5-4.  The  measure- 
ments of  Tisone  [Ti-73]  and  Barth  [Ba-66d]  are  in  qualitative  agree- 
ment with  discrepancies  of  less  than  a factor  of  two  over  a large  fraction 
of  the  common  altitude  range.  The  measurements  of  Meira  ]Me-71]  are 
distinctly  different  from  either  those  of  Tisone  [Ti-73  | or  Barth  |Ba-66d| 
by  the  pronounced  minimum  at  about  85  km.  Between  the  measurements 
of  Tisone  [Ti-73  | and  Meira  [Me-71 1 there  is  a difference  in  time  of 
day,  season  and  latitude:  the  measurements  of  Tisone  lTi-73|  were  ob- 
tained about  0500  LT  in  May  at  22  N whereas  those  of  Meira  [Me-71  | 
were  obtained  about  1400  LT,  in  January  and  February  at  38  N. 

There  are  four  sets  of  measurements  [AF-73,  RF-73ci,  TF-73,  ^ 
PB-74  I in  the  altitude  range  of  1 5 to  40  km;  these  are  plotted  on 
Fig.  5-4.  The  apparentlv  most  reliable  measurement  is  that  of  Patel 
et  al.  [PB-74|  in  which  the  measured  species,  NO,  is  unequivocallv 


67 


identified  by  use  of  a spin-flip  R:iman  laser.  This  measurement  clearly 
demonstrates  the  variation  of  nitric  oxide  density  between  sunrise  and 
noon  at  28  km.  The  most  extensive  set  of  measurements  are  those  of 
-Ackerman  et  al.  lAF-73)  based  on  absorption  by  NO  in  the  5.  2 |im  band 
usinjj;  the  setting;  sun  as  the  li^^ht  source.  As  the  solar  zenith  anc:le  in- 
creases above  90  , the  path  lietween  the  sun  and  the  detector,  at  40  km, 
passed  throuc;h  lower  and  lower  altitudes.  The  nitric  oxide  densities 
were  inferred  by  analysis  of  the  e.xtent  of  absorption  for  the  chane,in^ 
optical  paths.  The  measurements  of  Ridley  et  al.  |RS-73dj  depended 
on  the  reaction  of  ambient  nitric  oxide  with  ozone  in  the  detector  to  pro- 
duce the  chemiluminescence  which  was  measured.  Calibration  was  pro- 
vided by  an  on-board  source  of  nitric  oxide.  These  were  midday 
measurements  at  32  N latitude  in  March  at  altitudes  in  the  rane;e  of  17.4 
to  22.  9 km.  These  results  are  an  order  of  magnitude  smaller  than  the 
midday  results  of  Patel  et  al.  LPB-74]  at  28  km  in  the  same  latitude  re- 
gion in  October  of  the  same  year  (1973).  Measurements  by  Toth  et  al. 
lTF-731  provide  an  upper  limit  in  the  altitude  range  between  15  and  20 
km.  A ground- level  value  of  the  nitric  oxide  densitv,  recommended  tw 
Cadle  [Ca-74|,  is  plotted  on  Fig.  5-4. 

5.3.2  Model  Calculations 

The  model  calculations  of  nitric  oxide  densitv  in  the  lower  range 
of  altitudes  are  difficult  to  make  because  of  uncertainties  in  chemical 
and  transport  input  data.  As  in  the  case  of  atomic  oxygen,  a thorough 
comparison  of  the  various  models  was  not  made  for  this  review.  Sev- 
eral of  the  model  calculation  results  are  plotted  on  Fig.  5-4.  The  most 
e.xtensive  model  calculations  are  those  of  McConnell  and  McFlroy  [MM- 
731)  I and  Brasseur  and  Nicolet  lBN-73  1.  The  latter  investigators  have 
carried  out  a parametric  study:  the  results  which  are  shown  on  Fig.  5-4 
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are  for  calculations  using  the  upper  values  of  the  eddy  diffusion  coeffi- 
cients they  considered  and  the  variations  in  these  results  below  7 5 km 
corresptmd  to  nitrogen  atom  production  rates  of  0 (lower  cun’e)  and 
10  (upper  cuiwe)  atoms,  (cm^  sec).  Brasseur  and  Nicolet  [BN-73]  stress 
the  sensitivity  of  nitric  oxide  densities  to  values  of  the  eddy  diffusion 
coefficient.  For  a critique  of  eddy  diffusion  coefficient  values,  see 
Hunten  and  Strobel  [HS-74a|.  For  the  calculations  of  McConnell  and 
McElroy  [MM-73b]  two  models,  A and  D,  are  shown  on  Fig.  5-4. 

Model  A is  their  favored  model;  model  D differs  only  by  use  of  a lower 
quantum  yield  for  the  photolysis  of  HNO^.  The  value  of  the  quantum 
yield  is  presently  in  doubt  but  is  clearly  impoifant  for  the  calculation 
of  NO  densities,  particularly  in  the  altitude  regions  between  10  and 
30  km.  By  comparison  with  model  A,  model  D fits  the  obseiwational 
data  somewhat  better.  The  calculations  of  Hunt  [llu-73a]  and  Ciaitzen 
[Cr-7ll  appear  to  generally  overestimate  the  NO  densities.  The  calcu- 
lations of  Shimazaki  and  Laird  [SL-70a,  SL-72,  SL-72a]  at  0 and  60  N 
latitude,  bound  the  midlatitudes  of  interest  in  the  present  analysis.  The 
calculations  of  Shimazaki  and  Ogawa  [SO-74]  for  altitudes  below  50  km 
yield  densities  which  are  generally  not  in  accord  with  obseirations  and 
differ  by  factors  up  to  three  from  the  calculations  of  McConnell  and 
McElroy  [MM-73b].  The  earlier  calculations  of  Strobel  [St-72,  St-72a' 
are  shown  on  Fig.  5-4;  more  recent  calculations  using  an  improved  ver- 
sion of  Strobel's  calculations  [OJ-74]  are  shown  for  the  conditions  of  the 
obseiwations  of  Meira  [Me-71  | and  Tisone  [Ti-73].  (The  cuiwes  from 
these  calculations  [OJ-74]  as  drawn  on  F’ig.  5-4  are  accurate  onlv  at 
selected  altitudes  since  they  were  deduced  from  diurnal  nitric  oxide 
contours  given  [OJ-74]  only  for  a few  values  of  the  nitric  oxide  densitv.) 
There  is  qualitative  agreement  of  the  densitv  profiles  onlv  with  the  data 
from  Meira  ]Me-71  ].  In  the  altitude  range  between  75  and  100  km.  the 
calculated  values  are  smaller  than  anv  of  the  obseired  values  of  nitric 
oxide  density. 
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5.  3.  3 Selection  of  Nitric  Oxide  Profile 


fn  selecting  the  nitric  oxide  daytime  profile  below  120  km,  re- 
liance has  been  placed  exclusively  on  obseiwations.  Seasonal  and  lati- 
tudinal variations  are  not  taken  into  account  explicitly  although  these 
variations  may  be  contributing;  to  the  data  used.  In  the  region  below 
about  40  km,  nitric  oxide  is  no  longer  the  principal  odd  nitrogen  species: 
in  the  region  below  20  km,  NO2  and  HNO^  have  generally  greater  densi- 
ties than  nitric  oxide.  Use  of  the  selected  profile  in  calculations  which 
ignore  species  such  as  NO2  and  HNO^  will  have  little  meaning  at  these 
altitudes. 

In  the  70-115  km  inteiwal,  the  average  of  the  obseiwations  of 
Meira  [Me-7l],  Barth  LBa-66d]  and  Tisone  [Ti-73  ] were  used.  The 
observations  of  Pearce  [Pe-69a]  were  excluded  by  virtue  of  the  large 
discrepancy  between  these  obseiwations  and  others  for  this  altitude  re- 
gion as  well  as  between  60  and  70  km.  These  averages  are  given  in 
Table  5-2.  At  115  km,  the  density  was  obtained  bv  interpolation  using 
the  values  at  110  km  and  the  value  at  120  km  (see  Section  5.  2.  4). 

For  the  altitude  range  between  28  km  and  60  km,  the  datum  ol 
Patel  et  al.  lPB-74]  was  selected  as  the  absolute  NO  densitv  at  28  km: 
the  data  of  Ackerman  et  al.  [AF-73]  were  normalized  at  28  km  to  the 
datum  of  Patel  et  al.  [PB-74].  These  data  together  with  the  data  of  Hale 
[Ba-72a]  determine  the  nitric  oxide  profile  between  28  and  60  km.  See 
Table  5-2. 

Between  60  and  90  km,  a smooth  cuiwe  was  drawn  using  tlie 
previously  selected  values  at  60,  75  and  90  km  to  determine  the  locus 
of  the  cuiwe. 
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Table  5-2.  Selection  of  Values  of  Nitric  Oxide  Densities  Below  120  km. 


h. 

km 

INOI, 



cm 

<Me-71, 

Ba-66d 

Ti-73 

Kr-70a 

O 

Text* 

Ba-72a* 

R-S-73d 

115 

3.60(7) 

1 

110 

8. 1(7) 

3.  6(7) 

1.  5(7) 

4.35(7) 

4.  38(7) 

105 

9.6(7) 

4.6(7) 

2.75(7) 

5.65(7) 

100 

8.  2(7) 

5.  5(7) 

5.1(7) 

6.27(7) 

95 

6.0(7) 

6.0(7) 

5.  5(7) 

5.  83(7) 

90 

2.  8(7) 

6.0(7) 

5.  5(7) 

4.77(7) 

85 

1.35(7) 

6.0(7) 

5.  5(7) 

4.28(7) 

1 

i 

80 

1.  55(7) 

6.0(7) 

4.9(7) 

4.  15(7) 

1 

1 

75 

3.25(7) 

- 

3.4(7) 

3.33(7) 

70 

6.4(7) 

2.  2(7) 

4. 3(7) 

1 

1 

65 

5.  5(7) 

60 

1.4(8) 

55 

3.3(8) 

50 

5.  1(8) 

: 

45 

8.  5(8) 

40 

1.  25(9) 

35 

1.  75(9) 

30 

2.05(9) 

25 

1. 75(9) 

20 

1. 23(9) 

1. 75(8) 

15 

5.  8(8) 

10 

1.3(9) 

5 

3.4(9) 

0 

1.0(10) 

^See  text. 

Hale's  results  as  reported  in  Ba-72a. 
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At  25  km,  the  nitric  oxide  density  was  selected  Irom  the  ad- 
justed data  of  Ackerman  et  al.  |AF-73j;  at  20  km,  the  averaite  of  the 
datum  from  the  adjusted  cuiwe  of  Ackerman  et  al.  1AF-73|  and  of  that 
from  Ridley  et  al.  [Rt)-73dl  was  selected;  below  20  km,  a smooth  cuiwe 
was  drawn  using  the  selected,  20  km  datum,  a datum  for  16.  5 km  from 
the  adjusted  cuiwe  of  Ackerman  et  al.  lAF-73]  and  the  recommended 
[Ca-74]  nitric  oxide  density  at  ground  level.  Below  20  km,  where  only 
calculations  are  available  [MM-73b,  Le-72a.  Le-73),  peaks  in  the  cal- 
culated nitric  oxide  profile  jMM-73b]  have  been  ignored. 

The  nitric  oxide  profile  data  are  collated  in  Table  5-3  and 
shown  on  Fig.  5-5  for  altitudes  below  120  km. 

5.4  ADDITIONAL  DAYTIME  NITRIC  OXIDE  DATA 
5.  4.  1 Indirect  Determinations  of  Densities 

Nitric  oxide  densities  have  been  deduced  "indi rectiv" 

|Da-72e,  KB-74c,  NP-72a,  MR-74b,  01-74a,  ZS-74]  by  use  of  meas- 
urements on  electron  and  ion  densities  and  Lyman-alpha  radiation. 

Most  of  these  deductions  are  summarized  by  Mitra  and  Rowe  |MR-74b]: 

their  summary  demonstrates  that  between  65  and  100  km.  the  deduced 
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densities  of  nitric  oxide  have  values  between  10  and  10  cm  with  but 
few  exceptions.  This  summary  of  deduced,  nitric  oxide  densities  pro- 
vides [MR- 74b)  constraints  on  the  nitric  oxide  density  on  the  basis  of 
electron  densities  for  ambient  conditions,  for  conditions  resulting  from 
solar  activity  reversal,  for  flare  conditions  and  from  ion  composition 
measurements.  Not  included  in  this  summary  are  nitric  oxide  densities 
derived  from  rocket  experiment  data  |KB-74c,  ZS-74)  for  measuring 
ion  composition  and  Lvman-alpha  intensities  and  from  data  on  charged 
particle  density  variations  measured  during  an  eclipse  [01-74a).  For 
the  experiments  of  Zalpuri  and  Somavajulu  )ZS-74  ) between  85  and  90  km. 
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Table  5-3.  Nitric  Oxide  Profile  for  Noon  and 
Midnight  Conditions. 
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the  derived  densities  are  in  the  range  (1-4)  x 10  cm  whereas  belo\ 
85  km  there  is  a wide  disparity  in  deduced  nitric  oxide  densities  with 
the  noontime  profile  exhibiting  much  smaller  densities  than  the  early 
morning  profile. 

The  nitric  oxide  densities  deduced  by  Koshelev  and  Belinska 
[KB-74c]  lie  between  lo'^  and  2 x 10^  cm  ^ between  100  and  200  km. 
to  180  km,  these  deduced  nitric  oxide  densities  are  in  fair  agreement 
with  the  daytime  obseiwed  and  calculated  nitric  oxide  densities  shown 
Fig.  5-1,  provided  allowance  is  made  for  the  solar  activity.  Howeve 
above  180  km,  the  densities  of  Koshelev  and  Belinskaya  [KB-74c|  in- 
crease and  become  significantly  greater  than  expected  on  the  basis  ol 
the  treatment  of  Section  5.  2.  3.  The  densities  deduced  by  Oliver  [Ol- 
74a  1 between  110  and  280  km,  tend  to  be  slightly  lower  than  the  data 
calculational  results  shown  on  Fig.  5-1  for  altitudes  below  180  km,  1 
e.xhibit  a general  increase  above  180  km  in  contrast  to  direct  obseiwa 
tions  and  calculations. 

A difficulty  with  these  indirect  determinations  of  nitric  oxid( 
densities  is  the  necessity  of  modelling  correctly  the  production  and  1< 
mechanisms  for  nitric  oxide.  The  results  of  indirect  determinations 
nitric  oxide  densities  have  not  been  included  in  the  process  of  selecti 
a nitric  oxide  profile  in  the  present  repoid. 

Chakravai'ty  and  Mitra  [See  MR-74b]  have  deduced  nitric  ox 
densities  from  ion  composition  measurements  over  a solar  cycle  vai 
ation.  They  find  a corresponding  nitric  oxide  density  variation  of  a 
factor  of  four  down  to  100  km  with  a decreasingly  smalle r variation  ( 
to  90  km.  While  a smoothed  profile  (01-721  of  solar  flux  over  a sol; 
cycle  has  a larger  variation  in  Fj^  values  than  contained  in  the  da 
of  Table  5-1  above,  neveilheless,  the  results  of  Chakravai'ty  and  Mi 
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[MR-74b]  lend  credence  to  a variation  of  nitric  oxide  density  with  solar 
activity  below  120  km.  By  using  the  directly  obtained  data  at  110  km 
LBa-66d,  Me-71,  Kr-70a,  OJ-74,  Ti-73],  the  most  pi’obable  variation 
in  nitric  oxide  density  with  decimetric  solar  flux  is  found  to  be  a factor 
of  two  increase  between  values  of  90  and  145,  respectively:  the 

uncertainty  is  such  that  the  factor  could  lie  between  zero  and  four.  (Fur- 
ther discussion  is  presented  below  in  Section  5.  4.  2.)  However,  in  the 
present  const mction  of  profiles,  the  approximation  of  solar  flux  inde- 
pendence for  nitric  oxide  densities  below  120  km  is  retained  for  conven- 
ience and  simplicity  in  quantitative  description  of  the  profile. 

5.4,2  Late  Reports  on  Nitric  Oxide 

Two  reports  [GS-75a,  MS-75a]  have  recently  appeared  in  which 
nitric  oxide  densities  in  the  E- region  are  determined  indirectly  on  the 
basis  of  electron  and  ion  densities.  In  both  reports,  the  dependence  of 
the  deduced  nitric  oxide  density  on  solar  activitv  is  deivonstmted. 

Monro  and  Smith  [MS-75a]  have  found  at  130  km  that  the  nitric  oxide 
density  increases  with  the  sunspot  numi)er:  this  relation  is  shown  on 
Fig.  5-6  along  with  the  NO  densities  determined  directlv  or  calculated 
at  110  km  jBa-66d,  Me-71,  Kr-70a,  OJ-74,  Ti-73  ] and  which  have  been 
considered  in  previous  sections.  Thus  one  can  see  that  roughlv  the  same 
average  dependence  on  sunspot  numl)er  is  found  for  nitric  oxide  In’  direct 
measurement  or  calculation  at  110  km  as  bv  indirect  determination  at 
130  km  for  sunspot  numbers  greater  than  20.  The  absolute  concentra- 
tions of  nitric  oxide  determined  i)v  Monro  and  Smith  |MS-75a]  at  sun- 
rise, are,  however,  larger  by  factors  of  3 to  4 than  the  davtime  densitv 
measurements  of  Barth  |Ba-66c|  and  Meira  lMe-71  j.  These  results  als( 
disagree  with  calculations  [0.1-74,  Hu-73al.  The  method  of  analvsis 
used  by  Monro  and  Smith  [MS-75a|  requires  that  the  increase  in  electron 
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Fitr.  5-6.  Dpnpnrlpnrp  of  Nitrif  Oxirip  I')pnRitv  nn  Simsnnt  Nnmlmr 


density  in  the  sunrise  E-region  result  mainly  from  ionization  of  nitric 
oxide  by  Lyman-O!  radiation  for  a specified  range  of  solar  zenith  angles. 


Golshan  and  Sechrist  [GS-7  5a],  using  a different  set  of  ol)ser- 
vations  than  those  of  Monro  and  Smith  [MS-7 5a],  employ  the  obsen'ed 
[N0’^]/[02l  ratios  in  conjunction  with  atmospheric  composition,  flux 
and  cross-section  data  to  deduce  nitric  oxide  densities:  at  105  km,  they 
find  the  dependence  on  sunspot  number,  R^,  of  the  nitric  oxide  densitv 
as  shown  on  Fig.  5-6.  These  data  are  qualitatively  consistent  with  the 
data  at  110  km  shown  on  Fig.  5-6  and  are  I’oughly  a factor  of  two  larger. 
The  dependence  of  nitric  oxide  density  on  sunspot  number  is  related 
[GS-75a]  to  the  corresponding  variation  in  the  intensity  of  radiation  in 
the  10  to  130  A spectral  region.  Thus  the  rationale  for  variation  of  nitric 
oxide  density  with  solar  activity  differs  for  altitudes  near  100  km  and 
those  near  200  km.  From  the  data  assembled  on  Fig.  5-6,  a clear  jus- 
tification for  the  future  refinement  of  the  model  presented  in  this  review 
is  derivable. 

Golshan  and  Sechrist  [GS-7  5a]  also  find  a seasonal  dependence: 
above  100  km,  the  nitric  oxide  density  increases  from  winter  to  summer. 
This  is  in  contrast  to  the  calculations  of  Oian  et  al.  |OJ-74  | in  which 
the  calculations  for  winter  predict  more  nitric  oxide  than  those  for  sum- 
mer at  altitudes  above  100  km  and  for  a constant  molecular  oxygen 
density. 

A recent  calculation  on  the  seasonal  and  latitudinal  variation  of 
oxides  of  nitrogen  in  the  stratosphere  has  been  reported  [Vu-75|.  This 
work  consists  of  a preliminary,  two-dimensional  calculation  which  in- 
corporates photochemical,  veiiical  and  horizontal  (noidh-south)  transpoil 
processes.  The  preliminary  nature  is  primarily  a result  of  the  parame- 
terization of  the  ('zone  and  temperature  distributions  and  the  mean 
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meridional  motions  rather  than  their  calculations.  The  profiles  from 
this  calculation  at  45  latitude  for  both  summer  and  winter  are  shown  on 
Fig'.  5-4  for  altitudes  between  17  and  48  km. 

5.  5 NIGHTTI.MK  NITRIC  OXIDE  DENSITIES 
5.  5.  1 .Measurements 

Two  nighttime  rocket  flights  [OP-72,  OG-73]  have  gathered 
data  on  nitric  oxide  densities.  On  one  rocket  flight  (OG-73]  at  39  N 
latitude,  a ciwocooled  mass  spectrometer  was  the  detector;  on  the  com- 
panion flight  |OP-72l  an  uncooled  detector  was  used;  nevertheless  the 
results  of  the  two  measurements  [OP-72,  OG-73  [ agree  over  the  alti- 
tude range  between  170  and  210  km.  In  plotting  the  data  of  Offermann 
et  al.  [OP-72  I on  Fig.  5-1,  the  assumption  has  been  made  that  the  ratio 
of  the  M30  M28  signal  for  the  two  flights  [OP-72,  OG-73]  is  the  same 
as  the  ratio  of  the  densities;  MX  is  the  mass  of  species  X.  However, 
since  Offennann  and  Grossmann  [OG-73]  chose  not  to  consider  the  data 
obtained  with  the  uncooled  detector,  we  will  not  fuiiher  consider  the 
uncooled,  mass-spectrometer  data.  For  the  data  presented  by  Offer- 
mann and  Giossmann  [OG-73],  no  correction  was  made  for  the  possible 
contributions  to  the  nitric  oxide  signal  from  interactions  occurring  in 
the  mass  spectrometer;  thus  the  data  on  Fig.  5-1  represent  upper  limits. 
Nevertheless,  the  magiiitudes  of  the  deduced  densities  are  clearly  in 
contradiction  to  the  results  of  calculations,  as  discussed  below,  and  in- 
direct determinations  of  densities  [KB-74c]  since  the  nighttime  nitric 
o.xide  densities  of  Offermann  and  Grossmann  are  greater  than  those 
measured  or  calculated  for  daytime.  Since  these  data  cannot  be  ration- 
alized at  present,  the  selected  nighttime  density  profile  for  nitric  oxide 
will  relv  solely  on  calculations. 


81 


5.  b.  2 Calculations  and  Selection  of  Nitric  Oxide  Profile 


The  diurnal  variation  in  nitric  oxide  density  above  70  km  ob- 
tained from  the  calculations  of  Oran  et  al.  10J-74|  and  Hunt  [Hu-73a] 
shows  a smaller  density  for  nitric  oxide  at  midnight  than  at  noon  above 
an  altitude  of  100  km.  These  relative  values  are  also  obtained  in  the 
deduced  densities  of  Koshelev  and  Belinskaya  [KB-74cj.  In  the  calcu- 
lations [OJ-74,  Hu-73aJ,  the  difference  increases  with  altitude  and  at 
160  km,  the  noon  nitric  oxide  density  is  greater  than  the  midnight  den- 
sity by  factors  of  1.  8 [Hu-73a),  3.4  [OJ-74  for  Me-7l]  and  ~ 7.  1 

[OJ-74  for  Ti-73  |.  We  select  an  average  value  of  4 (= 

iNOl  . . . , J at  160  km,  a value  of  1 at  100  km  and  a maximum  value 
auidmght 

of  10  at  arbitrarily  high  altitudes;  then 


iNOl 


midnight 


INOl 

'^^°'Ml0(h  - 100)^  • 7200 


(0) 


where  h is  in  km. 

Below  70  km,  a significant  diurnal  variation  in  nitric  oxide 
density  is  predicted  by  Hunt  |Hu-73a|,  Shimazaki  and  Laird  (SL-70a, 
8L-72a|  and  Thomas  [Th-71al.  In  view  of  the  disparitv  between  these 
calculations  and  obseiwations  for  davdime  conditions  in  the  altitude 
range  below  70  km,  the  calculated  density  profiles  for  midnight  condi- 
tions must  be  viewed  skeptically.  Therefore,  the  fall  in  nitric  oxide 
densitv  below  70  km  will  have  a large  error  associated  with  it.  We  shall 
use  the  results  of  Shimazaki  and  Laird  |SL-70a,  SL-72|  from  70  to 
50  km  and  extrapolate  below  50  km.  The  use  of  these  calculations 
[SL-70a,  SL-72]  is  based  on  the  better  agreement  with  obseiwations  in 
the  region  between  50  and  70  km  than  foi'  the  calculations  of  Hunt  (Hu- 
73a  1.  The  values  are  shown  on  Table  5-3.  A partial  plotting  of  these 
values  is  given  on  Fig.  5-5.  Values  arc  interpolated  between  65  and 
85  km. 
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6.  DENSITY  PROFILE  OF  ATOMIC  NITROGEN 


6.  1 INTRODUCTION 

Despite  the  tact  that  a moderate  number  ol  measurements  and 
of  calculations  of  atomic  nitrogen  densities  have  been  made,  larger  dis- 
crepancies exist  among  the  measured  values  and  between  the  measured 
and  calculated  values  than  for  any  of  the  other  atmospheric  species  con- 
sidered in  this  report.  Thus  the  selected  profile  has  a large  unceidainty 
associated  uith  the  atomic  niti’ogen  densities. 

6.  2 MEASUREMENTS 

The  measurements  of  atomic  nitrogen  liave  been  made  at  alti- 
tudes above  100  km  on  rocket  flights  with  mass  spectrometers  [GH-68a. 
HN-72a,  MB-72b,  OP-72,  Po-72a]  and  by  a chemiluminescent  technique 
involving  the  reaction  of  nitrogen  atoms  with  ethylene.  The  data  trom 
these  measurements  are  shown  on  Fig.  6-1.  The  mass  spectromet ric 
measurements  are  difficult  to  make.  Important  corrections  to  the  raw 
data  must  be  made,  as  for  example,  to  account  for  the  ionization  of 
molecular  nitrogen  to  N’2  and  N [HN-72a].  Also,  the  possible  influ- 
ence of  vibrationally  excited  nitrogen  has  been  considered  [DB-74]. 

One  may  also  ask  whether  electron  excitation  of  molecular  nitrogen  to 
predissociating  states  contributes  ultimately  to  the  measured  nitrogen 
ion  signal.  At  electron  energies  of  25  eV,  for  example,  the  cross 
section  for  dissociation  of  molecular  nitrogen  l\\'i-66|  (accounted  for 
primarily  iiy  i)redissociation  of  excited  states  of  molecular  nitrogen)  is 
comparable  to  the  cross  section  for  ionization  of  molecular  nitrogen. 

Can  the  nitrogen  atoms  formed  by  electron  impact  contriliute  to  the  nitro- 
gen ion  signal  and  in  significantlv  different  wavs  in  the  atmosphere  and 
the  laboratoiT  caliiiration  measurements? 
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The  chemiluminescent  measurements  were  based  on  a rate  co- 
efficient derived  from  experiments  in  one  laboratorv  [BB-74|.  On  the 
^ other  hand,  any  error  in  the  rate  coefficient  is  more  likely  to  result  in 

too  large  a rate  coefficient  which  would  have  the  effect  of  increasing  the 
nitrogen  atom  density  calculated  from  the  proper  coefficient. 

’ The  measurements  shown  on  Fig.  6-1  were  conducted  in  the 

early  morning  or  at  night  with  one  exception.  There  is  no  significant 
correlation  of  the  spread  in  measurements  with  latitude;  the  experiments 
conducted  at  Volgograd  (48  N)  have  yielded  larger  nitrogen  atom  densi- 
ties tlian  e.xperiments  conducted  elsewhere.  Note  also,  that  the  e.xperi- 
meats  of  Hickman  and  Nier  |HN-72aj  yield  only  upper  limit  estimates 
of  atomic  nitrogen  densities.  For  the  experiments  of  Hickman  and  Nier 
|HN-72a|  and  Offermann  et  al.  [OP-721,  we  have  used  the  CIR.V  1972 
[CI-72]  mean  atmosphere  to  convert  from  the  signal  ratios  of  mass  14  to 
mass  28  to  nitrogen  atom  densities  and  for  the  experiment  of  Pokhunkov 
r lPo-72a|,  a mean  of  the  atmospheres  corresponding  to  T^  - 1200  K and 

1300  K has  been  used  for  the  conversion.  The  conversion  of  the  data  of 
Martvnkevich  and  Bvuro  |MB-721)1  has  been  given  bv  Balabanova  et  al. 

, |BB-74l. 

The  measurements  shown  in  I'ig.  6-1  have  been  examined  with 
respect  to  the  solar  10.  7 cm  flux  observed  for  the  day  of  the  rocket 
; flight:  any  correlation  of  the  atomic  nitrogen  density  with  the  10.7  cm 

flux  is  obscured  by  differences  in  the  atomic  nitrogen  densitv  of  one  to 
j two  orders  of  magnitude  for  measurements  corresponding  to  verv  simi- 

lar solar  flux  values. 
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6.3  MODEL  CALCULATIONS 

The  results  of  model  calculations  for  atomic  nitrogen  are 
shown  on  Fig.  6-2.  The  calculations  of  Hunt  lHu-73a|  and  Strobel  and 
colleagues  [St-71,  St-72a,  OJ-74]  have  included  ionic  reactions  as  well 
as  neutral  reactions  and  transport.  The  calculations  of  Shimazaki  and 
Laird  [SL-70a,  SL-72,  SL-72a|  have  not  included  ionic  reactions  and 
their  calculations  of  atomic  nitrogen  densities  are  in  error  at  least  for 
altitudes  above  85  km.  Nicolet  |Ni-70]  has  calculated  atomic  nitrogen 
densities  but  with  photo-equilibrium  conditions.  The  calculations  of 
Hunt  [Hu-73a|  are  more  comparable  with  the  Oran  et  al.  (OJ-74  | cal- 
culations for  the  e.xperimental  conditions  of  Tisone  [Ti-73]  than  for 
those  of  Meira  [Me-71]  on  the  basis  of  solar  activity;  the  Hunt  calcula- 
tions yield  densities  at  least  a factor  of  three  smaller  than  those  of 
Oran  et  al.  for  conditions  of  Tisone's  e.xperiment.  The  calculations  of 
Oran  et  al.  [OJ-74  ( show,  for  both  day  and  night  conditions,  a larger 
nitrogen  atom  density  for  the  conditions  of  Tisone’s  experiment  (Ti-73| 
than  for  the  condition  of  Meira's  e.xperiment  (Me-71]  at  altitudes  above 
100  km.  This  is  consistent  with  the  expectation,  based  on  the  analysis 
of  Section  5.  2.  3,  that  a lower  solar  flux  would  result  in  a larger  atomic 
nitrogen  density.  However,  for  the  higher  altitudes,  the  effect  of  solar 
flux  on  density  is  apparently  considerat)ly  smaller  for  atomic  nitrogen 
tlian  for  nitric  oxide,  at  least  in  the  daytime. 

6.4  SELECTION  OF  ATOMIC  NITROGEN  PROFILE 

Because  of  the  verv  large  disparity  among  measurements  and 
calculations,  a noon  and  midnight  profile  will  be  deiived  by  freely  draw- 
ing a cuiTC  through  the  calculated  profiles  and  no  attempt  will  be  made 
to  consider  the  effect  of  solar  activity  on  atomic  nitrogen  densities.  The 
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greatest  weight  is  given  to  the  recent  calculations  of  Oran  et  al.  lOJ-74]: 
zero  weight  is  given  to  the  measurements  although  this  might  be  negligent. 

Data  for  the  selected  profile  are  listed  in  Table  6-1  and  shown 
on  Fig.  6-3  for  altitudes  up  to  230  km. 
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Above  230  km,  [N  | (jN],  l^'2l^230  ^^2^  tke  subscript  ind 

the  nitroe;cn  atom  to  nitrogen  molecule  ratio  at  230  km. 


Table  6-1. 

Nitro<]:en  Atom  Densities  for  Noon  and 
Alidnight  Conditions. 
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7.  DENSITY  PROFILE  OF  NITROGEN  DIOXIDE 


7.  1 INTRODUCTION 

Only  recently  has  nitrogen  dioxide  been  detected  at  al 
much  above  ground  level  [MK-68b,  \IM-69,  GM-70c].  Thus  I 
is  not  surprising  that  nitrogen  dioxide  densities  show  a scattei 
an  order  of  magnitude  in  the  principal  altitude  region  of  the  mi 
ments,  i.  e.  , the  stratosphere. 

The  selected  nitrogen  dioxide  profile  applies  to  the  al 
region  from  ground  level  to  160  km.  This  profile  is  based  on 
ments  and  calculations  for  altitudes  below  30  km  but  above  thi; 
model  calculations  are  the  sole  basis. 

7.2  DAYTIME  NITROGEN  DIOXIDE  DENSITIES 
7.  2.  1 Measurements 

The  measurements  of  nitrogen  dioxide  are  generallv 
to  the  altitude  range  from  about  14  to  34  km  and  in  that  range 
large  scatter.  The  scatter  in  the  measured  densities  is  paidii 
large  only  in  the  range  from  about  14  to  18  km.  This  is  showi 
P’ig.  7-1.  Note  that  the  densities  attributed  to  Farmer  [Fa-7‘ 
based  on  a tentative  value  of  the  mixing  ratio  determined  for  t 
15  to  20  km.  By  using  the  extensive  data  on  nitric  acid.  Schif 
has  calculated  a mixing  ratio  band  (converted  here  t('  densitie: 
the  U.  S.  Standard  .Atmosphere  [US-62  1 number  densities)  for 
tude  range  from  14  to  27.  6 km;  the  derived  band  of  nitrogen  d 
densities  favors  small  values  Ijelow  20  km.  .Also  supporting  t 
nitrogen  dioxide  densities  are  the  values  derived  by  Schiff  [Sc 
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the  data  of  Ackerman  et  al.  lAF-74  1,  the  data  of  Fontanella  et  al.  (FG 
74]  as  reported  by  Hard  [Ha-74],  and  the  value  calculated  bv  Murcra; 
based  on  data  of  Ackerman  et  al.  [AF-73a]  and  as  reported  by  Hard 
[Ha- 7 4]. 

The  measurements  of  nitrogen  dioxide  densities  have  general 
been  absorption  measurements  with  the  sun  as  the  radiation  source, 
the  measurements  shown  on  Fig.  7-1  but  with  the  three  exceptions  not 
below  are  from  observations  of  the  N02"^'2  6.  25  pm:  this  ban 

is  strong  and  relatively  free  from  interference.  Heretofore,  a difficc 
with  the  use  of  the  '^^.nd  has  been  uncei'tainty  about  the  rotational 
quantum  assignments.  Reference  to  new  laboratory  data  of  Mure ray  < 
the  1^2  band  of  nitrogen  dioxide  is  made  by  Schiff  [Sc-74bl;  Schiff  has 
reported  the  density  at  30  km  based  on  these  new  data  as  calculated  Ir 
Murcray  and  as  shown  on  Fig.  7-1.  The  measurements  of  .A.ckerman 
and  Muller  [AM-72,  AM-73]  and  of  Farmer  [Fa-74a]  have  used  the 

-r'o  band  at  3.4  pm  which  is  much  weaker  than  the  band  and  is 
1 o 

overlapped  by  methane  bands.  Harries  [Ha-73e]  has  used  high  resol 
emission  measurements  in  the  sub-millimeter  infrared  spectral  regi( 
The  datum  of  Harries,  as  given  on  Fig.  7-1,  is  a revision  of  earlier 
nitrogen  dioxide  densities  [Ha-73e];  the  revised  value  is  leported  by 
Farmer  jFa-74a]  as  well  as  Schiff  lSc-74b]. 

AU  of  the  measurements  of  Fig.  7-1  luave  been  made  l>etwee 
about  30  N and  50  N latitude,  except  those  of  Harries  [Ha-73e]:  for 
these  aircraft  obseiwations,  the  Latitude  and  longitude  varied  widelv  1 
40  S to  48  N and  from  8 W to  151  E.  The  accuraev  of  the  measure:-: 
was  limited  and  anv  dependence  on  geographic  location  is  obsciii  o-i 
the  measurements  of  Fig.  7-1  were  made  for  zenith  ancles  cl 
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Additional  measurements  not  given  on  Fig.  7-1  have  been  made 
by  Brewer  et  al.  [BM-73a];  these  measurements  have  not  been  included 
on  the  basis  of  the  discussion  given  by  Schiff  (Sc-74b], 

7.2.2  Model  Calculations 

The  model  calculation  results  are  shown  both  on  Fig.  7-1  and 
on  Fig.  7-2.  For  the  calculations  of  McConnell  and  McElroy  [MM-73b), 
the  cuiwes  on  Figs.  7-1  and  7-2  represent  their  model  A which  is  based 
on  the  best  estimate  of  jDarameters;  the  calculations  apply  to  30  N lati- 
tude at  equinox.  The  calculations  of  Shimazaki  and  Ogawa  lSO-74c], 
which  are  shown  on  Figs.  7-1  and  7-2,  are  the  model  4 calculations  of 
these  investigators  (see  Table  5 of  SO-74c)  and  apply  to  45  N latitude  at 
equinox.  The  shape  of  the  profile  is  similar  to  that  of  McConnell  and 
McElroy  lMM-73bJ  but  the  profile  has  smaller  densities;  these  dif- 
ferences may  be  due  to  differences  in  rate  coefficients  lSO-74cl.  A 
third  calculation,  applicable  at  low  altitudes,  was  carried  out  by  Crutzen 
[Cr-74a]  for  sec  6 = 1.4  where  0 is  the  solar  zenith  angle.  This 
calculation  is  in  best  agreement  with  experimental  data  between  25  and 
30  km  and  also  in  good  agreement  with  the  data  of  Murcray  lHa-74]  at 
16  km.  The  calculations  of  Hunt  lHu-73a]  yield  densities  an  order  of 
magnitude  larger  at  60  km  than  those  calculated  by  McConnell  and 
McElroy  [MM-73b].  The  calculations  of  Shimazaki  and  Laird  [SL-72a] 
have  larger  nitrogen  dioxide  densities  between  40  and  60  km  than  other 
calculations  and  smaller  densities  between  60  and  70  km. 

7.2.3  Selection  of  Nitrogen  Dioxide  Profile 

The  pi-ofile  is  selected  mainly  on  the  basis  of  calculations  but 
taking  into  account  the  experimental  data.  At  ground  level,  the  datum 
reported  by  Cadle  lCa-74l  is  used.  Between  ground  level  and  10  km, 
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the  calculations  of  McConnell  and  McElroy  [MM-73b]  are  adopted;  from 
10  to  near  40  km,  the  calculations  of  Cnjtzen  [Cr-74a]  are  used  and 
from  above  40  to  55  km,  those  of  McConnell  and  McElroy  [MM-73b]  and 
for  60  km  and  above,  those  of  Hunt  [Hu-73a].  Interpolation  between  the 
results  of  Crutzen  [Cr-74a]  at  35  km  and  those  of  McConnell  and 
McElroy  [MM-73b]  at  4 5 km  and  between  the  results  of  McConnell  and 
McElroy  [MM-73b]  at  50  km  and  those  of  Hunt  [Hu-73a|  at  65  km  pro- 
vided the  transitions  from  one  calculation  to  the  other.  Above  70  km,  the 
profile  of  Hunt  [Hu-73a]  has  been  smoothed.  Table  7-1  and  Fip;.  7-3 
show  the  selected  data. 

7.2.4  Late  Reports  on  Nitrogen  Dioxide 

After  the  nitrogen  dioxide  profile  had  been  selected,  a calcula- 
tion on  the  seasonal  and  latitudinal  variation  of  oxides  of  nitrogen  in  the 
stratosphere  was  reported  [Vu-75]  as  indicated  in  Section  5.4.  2 above. 
The  profiles  from  this  calculation  for  45°  latitude  both  summer  and 
winter  are  shown  on  Figs.  7-1  and  7-2.  Between  41  and  60  km,  these 
profiles  [Vu-75]  are  in  agreement  with  other  profiles  shown  on  Figs. 

7-1  and  7-2  but  below  41  km  show  generally  smaller  nitrogen  dioxide 
densities.  The  calculation  of  Vupputuri  [Vu-75]  exhibits  a decline  in 
the  nitrogen  dioxide  mixing  ratio  from  the  summer  to  the  winter  pole  for 
altitudes  below  about  41  km  and  this  is  indicated  by  the  curves  on  Figs. 
7-1  and  7-2  (densities  were  computed  from  the  calculated  mixing  ratios 
by  using  the  U.  S.  Standard  Atmosphere  [US-66]).  A superficial  assess- 
ment of  the  measurements  shown  on  Fig.  7-1  does  not  support  the  sea- 
sonal dependence  of  nitrogen  dioxide  densities  as  calculated  by  Vupputuri 
[Vu-75];  however  both  the  calculations  and  the  measurements  are  in  the 
early  stages  of  development. 
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Table  7-1.  Nitrogen  Dioxide  Densities  for  Noon  and  Midnight  Conditions. 


h, 

km 

[NO^],  cm”^ 

m 

[NOg],  cm"^ 

Noon 

Midnight 

Noon 

Midnight 

0 

2.  5(10) 

3.  5(10) 

85 

1.2(4) 

1.2(4) 

5 

8.3(9) 

1. 2(10) 

90 

6.4(3) 

6.4(3) 

10 

1.4(9) 

2.7(9) 

95 

3.4(3) 

3.4(3) 

15 

1.4(9) 

2.0(9) 

1.8(3) 

1.8(3) 

20 

1.  8(9) 

2.  5(9) 

1.1(3) 

1.1(3) 

25 

2.4(9) 

4.15(9) 

no 

6.7(2) 

6.  7(2) 

30 

2.  5(9) 

4.  55(9) 

115 

4.3(2) 

4.3(2) 

35 

1.  25(9) 

3.0(9) 

120 

2.  8(2) 

2.  8(2) 

40 

3.4(8) 

1.6(9) 

125 

1.9(2) 

1.9(2) 

45 

7.1(7) 

9.2(8) 

130 

1.4(2) 

1.4(2) 

50 

7.8(6) 

5.  2(8) 

135 

1. 15(2) 

1.  15(2) 

55 

2.3(6) 

3.0(8) 

140 

9.  5(1) 

9.  5(1) 

60 

7.0(5) 

1.4(8) 

145 

8.  0(1) 

8.  0(1) 

65 

2.6(5) 

5.  5(7) 

150 

7.0(1) 

7.0(1) 

70 

1.0(5) 

1.2(7) 

155 

6.0(1) 

6.0(1) 

75 

5.0(4) 

1.0(6) 

160 

4.6(1) 

4.6(1) 

80 

2.4(4) 

3.0(4) 
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7.3  NIGHTTIME  NITROGEN  DIOXIDE  DENSITIES 
7.3.1  Model  Calculations 

Only  model  calculations  are  available  for  estimating  the  night- 
time nitrogen  dioxide  densities. 

The  calculations  of  Shimazaki  and  Ogawa  [SO-74c]  show  the 
diurnal  variation  in  nitrogen  dioxide  density  at  20  km.  Much  of  the 
nitric  oxide  is  converted  to  nitrogen  dioxide  within  a short  time  inteiwal 
about  twilight  and  the  nitrogen  dioxide  density  is  constant  until  sunrise. 
Wliile  other  species  densities  increase  at  sunset  and  during  the  night, 
i.  e.  , NOg  and  NgOg,  these  densities  are  small.  The  calculated  ratio 
of  midnight  to  noon  NO2  densities  is  1.  87;  in  combination  with  the 
selected  nitrogen  dioxide  daytime  profile,  this  ratio  would  result  in  a 
nitrogen  dioxide  nighttime  density  which  exceeds  the  total  of  the  daytime 
densities  of  NO  and  NO2  from  the  corresponding  selected  profiles.  This 
difficulty  apparently  reflects,  in  paid,  the  relatively  small  values  of  the 
daytime  nitrogen  dioxide  density  calculated  by  Shimazaki  and  Ogawa 
[SO-74cJ  at  20  km  in  comparison  to  the  measured  values. 

Below  20  km,  the  sum  of  the  densities  of  nitric  oxide  and  ni- 
trogen dioxide  corresponds  to  a fixed  value  of  the  mixing  ratio  [Le-73]: 
at  night,  nitric  oxide  is  rapidly  converted  to  nitrogen  dioxide  lLe-73]. 

At  60  km,  the  calculations  of  Hunt  [Hu-73a|  and  of  Shimazaki 
and  Laird  [SL-70a,  SL-72,  SL-72a]  show  the  nighttime  nitrogen  dioxide 
density  to  be  equal  to  the  daytime  nitric  oxide  density  (plus  the  daytime 
nitrogen  dioxide  density,  presumably)  and  this  obtains  also  between  40 
and  60  km  according  to  the  calculations  of  Shimazaki  and  Laird  [SL-70a, 
SL-72,  SL-72a].  At  60  km,  the  rati'''  of  midnight  to  noon  nitrogen  diox- 
ide densities  given  by  Hunt  (Hu-73a]  is  794,  by  Shimazaki  and  Laird 
[SL-70a,  SL-72,  SL-72a]  is  1232  and,  by  use  of  the  nitric  oxide  and 


99 


r 


nitrogen  dioxide  daytime  profiles  of  Figs.  5-5  and  7-3,  is  212.  For  the 
case  of  the  Hunt  calculation  [Hu-73a],  the  large  ratio  results  primarily 
from  a greater  calculated  daytime  nitric  oxide  density  at  60  km  than 
indicated  by  observations  on  which  the  profile  of  Fig.  5-5  is  based.  In 
the  case  of  the  Shimazaki  and  Laird  calculation  [SL-70a,  SL-72,  SL- 
72a],  the  large  ratio  results  primarily  from  a smaller  calculated  day- 
time nitrogen  dioxide  density  at  60  km.  Differences  in  calculated 
daytime  nitric  oxide  and  nitrogen  dioxide  densities  at  high  latitudes 
[SL-72a]  also  result  in  larger  ratios  of  midnight  to  noon  nitrogen  dioxide 
densities  than  212  as  found  from  Figs.  5-5  and  7-3. 

Above  80  km,  the  calculations  of  Hunt  [Hu-73a]  show  little 
diurnal  variation  in  nitrogen  dioxide  densities.  This  is  expected  on  the 
basis  that  reactions  of  atomic  oxygen  dominate  the  formation  and  loss  of 
nitrogen  dioxide  above  80  km.  The  absence  of  a diurnal  variation  is  also 
indicated  by  the  calculations  of  Shimazaki  and  Laird  [SL-72a|  whereas 
in  earlier  calculations  for  lower  latitudes  [SL-70a,  SL-72]  the  diurnal 
variation  persists  to  90  km. 

7.3.2  Selection  of  Nitrogen  Dioxide  Profile 

The  selected  midnight  nitrogen  dioxide  profile  is  presented  in 
Table  7-1  and  Fig.  7-3.  Between  ground  level  and  60  km,  the  daytime 
nitric  oxide  is  converted  entirely  to  nitrogen  dioxide  and  the  sum  of  this 
quantity  and  the  daytime  nitrogen  dioxide  density  is  taken  to  be  the  mid- 
night nitrogen  dioxide  density.  This  results  in  a ratio  of  midnight  to 
noon  nitrogen  dioxide  densities  of  1.  4 at  20  km  as  compared  to  the  value 
of  1.  87  as  computed  by  Shimazaki  and  Ogawa  [SO-74c].  Above  80  km, 
any  diurnal  variation  in  NO2  is  neglected.  Between  60  and  80  km, 
nitrogen  dioxide  midnight  densities  are  obtained  by  interpolation  guided 
by  the  relative  change  in  the  ratio  of  midnight  to  noon  nitrogen  dioxide 
density  values  calculated  by  Hunt  [Hu-73a]  in  this  altitude  range. 
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8.  CARBON  DIOXIDE  PROFILES 


8.  1 INTRODUCTION 

Carbon  dioxide  is  apparently  well  mixed  in  the  atmosphere  up 
to  altitudes  near  100  km.  There  is  extensive  evidence  for  this  up  through 
the  troposphere  [Ca-74],  at  altitudes  in  the  range  44-62  km  [SH-69a] 
and,  by  deduction  from  measurements  of  CO2  densities  in  the  lower 
thermosphere  [OG-73],  up  to  100  km.  Above  100  km,  CO2  is  generally 
taken  to  be  in  diffusive  equilibrium  on  the  basis  of  observations  [OG-73, 
OP-72,  OZ-71]. 

8. 2 MEASUREMENTS 

Measurement  of  CO2  densities  up  to  the  stratosphere  are  sum- 
marized by  Cadle  [Ca-74].  The  annual  mean  mixing  ratio  is  reported 
to  be  3.22  x 10”^  [Ca-74]  for  a 45°N  latitude.  Between  11  and  20  km, 

_4 

the  value  recommended  is  3.  21  x 10  . Between  20  and  37  km,  Acker- 

-4 

man  [Ac-74]  found  the  volume  mixing  ratio  to  be  2.  6 ±0.  7 x 10  by 
absorptiometry:  between  44  and  62  km,  Scholz  et  al.  [SH-69a]  found 

_4 

the  mixing  ratio  to  be  (3.  228  ±0.  025)  10  . Observations  in  the  lower 

thermosphere,  at  altitudes  between  120  and  140  km  have  been  made  on 
rocket  flights  [OG-73,  OP-72,  OZ-71]  using  a helium-cooled  mass 
spectrometer.  From  these  measurements,  the  densities  of  CO2  at 
three  altitudes,  124,  129,  and  135  km,  have  been  obtained.  These 
measurements  are  consistent  with  the  assumption  of  diffusive  equilibri- 
um above  100  km.  Recent  measurements  [SU-74  | of  the  infrared  emis- 
sion from  the  band  of  CO2  should  yield  a CO2  profile  between  60  and 
150  km  when  reduction  of  the  data  is  completed. 


3 
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All  the  observations  considered  here  are  plotted  on  Fi<^.  8-1 
in  terms  of  the  volume  mixing  ratios.  Also,  the  rocket  flight  measure- 
ments are  plotted  as  CO2  densities. 

8.  3 MODEL  CALCULATIONS 

Three  model  calculations  for  CO2  have  been  recently  made 
[Hu-73a,  WlVI-72,  HO-70].  The  most  extensive  of  these  calculations  is 
that  of  Hunt  [Hu-73a].  The  model  considered  an  oxygen- hydrogen- 
nitrogen- carbon  atmosphere  between  60  and  160  km  with  both  molecular 
and  eddy  diffusion  for  the  long-lived  species  including  CO2.  Photolvtic 
and  chemical  reactions  for  CO2  were  included.  Onlv  at  the  upper 
boundary  was  diffusive  equilibrium  assumed;  thus  it  is  worth  noting  that 
the  calculated  CO2  densities,  at  least  above  120  km,  are  essentiallv 
identical  to  the  densities  calculated  on  the  basis  of  diffusive  equilibrium 
above  120  km  but  starting  with  the  calculated  CO^  density  lHu-73a]  at 
120  km. 

The  calculations  of  Wofsy  et  al.  [WM-72|.  which  e.xtend  only 
to  about  120  km,  agree  with  those  of  Hunt  [Hu- 73a]  (see  mixing  ratio  plot 
of  Fig.  8-1)  except  below  about  100  km  where  slightly  different  values 
of  the  mixing  ratio  were  used  in  the  two  investigations. 

The  earlier  calculations  of  Hays  and  Olivero  [HO-70|  are  also 
shown  on  Fig.  8-1  on  both  mixing  ratio  and  density  profile  plots.  The 
densities  have  been  calculated  from  the  mixing  ratio  by  using  the  CIIIA- 
65  mean  atmosphere  model  [CI-65].  These  calculations  yield  smaller 
values  than  the  calculations  of  Hunt  [Hu-73a]  for  the  mixing  ratios  and 
CO2  densities.  The  chemical  system  of  Hays  and  Olivero  [HO-70]  was 
restrictive  bv  comparison  with  that  of  Hunt  [Hu-73a|  and  according  to 
Hunt  [Hu-73a],  Hays  and  Olivero  [HO-70]  specified  O and  OH  profiles 
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Fig.  8-1.  Carbon  Dioxide  Density  and  Mixing  Ratio  Profiles 
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rather  than  calculated  them.  In  addition  to  the  smaller  calculated  CO2 
densities,  the  Hays  and  Olivero  [HO-70]  calculations  predict  a departure 
from  mixing  above  80  km  in  contrast  to  the  value  of  100  km  for  the  cal- 
culations of  Hunt  [Hu-73aJ  and  Wofsy  et  al.  [WM-72]. 

The  calculations  of  Hunt  [Hu-73a]  were  made  for  diurnal  con- 
ditions. Only  a small  diurnal  variation  was  reported;  thus  at  100  km, 
CO2  increased  by  25%  in  daytime  over  the  nighttime  value. 

8.4  SELECTION  OF  CO2  PROFILE 

By  comparing  the  model  calculations  and  measurements  in  the 
density  plot  of  Fig.  8-1  for  altitudes  above  100  km,  one  finds  the  meas- 
ured CO2  densities  to  be  significantly  smaller  than  the  calculated  densi- 
ties. Both  the  measurements  and  calculations  indicate  diffusive 
equilibrium  above  100  km.  The  measurements,  however,  are  limited 
to  three  specific  altitudes,  at  present.  We  choose  to  give  equal  weight 
to  the  measurements  [OG-73)  and  to  the  model  calculation  of  Hunt 

[Hu-73a]  for  120  km  and  higher  altitudes.  Thus  we  select  a CO„  density 
7-3 

of  5.  65  X 10  cm  at  120  km  and  calculate  densities  at  higher  altitudes 
by  assuming  diffusive  equilibrium.  At  altitudes  of  100  km  and  below, 
complete  mixing  is  assumed  in  accordance  with  the  two  recent  model 
calculations  [Hu-73a,  WM-72)  and  observations  below  63  km  [SH-69a, 
Ca-74];  a value  of  3.  2 x 10  is  recommended  for  the  mixing  ratio.  In 
the  interval  between  100  and  120  km,  values  of  the  CO2  density  are  ob- 
tained by  interpolation  as  shown  on  Fig.  8-2.  Also  shown  on  Fig.  8-2 
are  the  densities  calculated  by  Hunt  lHu-73al,  by  use  of  the  ATMOSU 
program  [HS-75],  and  data  of  Offermann  and  Grossmann  [OG-73  ].  The 
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ATMOSU  program  lHS-75|  uses  a mixing  ratio  of  3.2  X 10  . The 

recommended  mixing  latios  and  COg  densities  are  summarized  in  Ta- 
ble 8-1.  Any  diurnal  variation  is  neglected  in  the  recommended  profile. 

8.  5 COMPARISON  OF  PROFILES 

The  profile  recommended  above  differs  at  altitudes  of  120  km 
and  greater  from  a recently  proposed  profile  [SD-73b|.  This  difference 
arises  as  a result  of  the  manner  of  selecting  the  density  of  CO2  at 
120  km  from  the  data  of  Offermann  and  coworkers.  In  the  reports  of 
the  German  group  [OG-73,  OP-72]  the  ratios  of  the  CO2  density  to  the 
(measured)  Ar  density  were  given.  The  CO2  density  can  be  obtained 
from  these  ratios  by  using  the  measured  Ar  density  as  reported  in  Offer- 
mann and  Grossmann  [OG-73  |;  this  was  the  procedure  used  above. 
However,  for  the  profile  of  Stephens  et  al.  [SD-73b],  the  CO2  density 
was  obtained  by  multiplying  a value  from  the  curve  representing  the  re- 
ported density  ratios  [OP-72]  by  the  density  of  Ar  as  given  in  CIRA-65; 


Table  8-1.  Carbon  Dioxide  Mixing  Ratio  and  Densities. 


h,  km 

Mixing  Ratio 

_3* 

[CO2],  cm 

0 to  100 

3.  2 X 10“^ 

105 

1.30(9) 

no 

4.8(8) 

115 

1.7(8) 

120 

5.  65(7) 

> 120 

Assume  diffusive  equilibrium 

*No  diurnal  variation;  values  to  be  used  for  noon  and  midnight. 
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this  procedure  yields  a CO„  density  at  120  km  4.  5 times  larger  than  the  i 

CO„  density  as  obtained  above  and  shown  in  Fig.  8-1  (i.  e.  , the  Ar  den- 
sity  in  ClRA-65  is  about  4.  5 times  the  measured  argon  density  [OG-73]).  j 

As  a consequence  of  these  considerations,  the  CO2  density  used  in  the  ^ 

recommended  profile  of  Stephens  et  al.  [SD-73b]  is  too  large  provided 
the  measured  COg  densities  [OG-73,  OP-72]  are  accepted  as  being 
accurate.  The  actual  difference  at  120  km  in  the  COg  profile  recom- 
mended in  the  present  report  and  that  recommended  by  Stephens  et  al. 

[SD-73b]  is  smaller  than  corresponds  to  the  factor  of  4.  5 as  a result  of 
the  averaging  procedure  used  here  to  obtain  the  recommended  density 
of  CO2  at  120  km;  the  ratio  of  the  COg  density  at  120  km  in  the  profile 
of  Stephens  et  al.  [SD-73b]  to  that  used  here  is  only  1.  9 rather  than  4.  5. 
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9.  DENSITY  PROFILE  FOR  WATER  VAPOR 


9. 1 INTRODUCTION 

Water  density  measurements  and  calculations  are  presented 
for  altitudes  up  to  160  km.  The  experimental  measurements  are  limited 
to  altitudes  below  about  60  km.  For  the  data  of  the  experiments,  there 
is  a considerable  variability.  The  variability  is  considered  by  some 
[Ca-74]  to  merit  further  consideration  while  others  [HS-74a]  dismiss 
the  higher  set  of  water  density  values  as  resulting  from  water  contami- 
nation or  other  errors  in  the  experiments.  Continued  improvements  in 
the  experimental  techniques,  particularly  those  based  on  emission  or 
absorption  spectroscopy  (which  has  only  recently  been  extensively  em- 
ployed), are  most  likely  to  resolve  the  different  appraisals. 

9. 2 MEASUREMENTS 

A discussion  of  measurements  has  been  given  by  Cadle  [Ca-74] 
for  water  vapor  at  the  surface,  in  the  troposphere  and  in  the  stratosphere. 
A profile  of  water  vapor  mixing  ratios  from  ground  to  16  km  prepared 
by  Sissenwine  on  the  basis  of  radiosonde  data  [GS-66b|  is  reported  by 
Cadle  [Ca-74]  for  nominal,  midlatitude  conditions.  This  profile  is 
shown  on  Fig.  9-1.  For  comparison,  measurements  of  Sissenwine  et 
al.  |SG-68a,  SG-68b]  and  an  older,  midlatitude  mean  profile  by  Gutnick 
[Gu-61j  are  also  shown  on  Fig.  9-1  over  the  common  altitude  range. 

Note  that  values  of  the  mixing  ratio  used  here  are  in  parts  per  million 
by  volume  (ppmv)  whereas  in  the  literature  for  water  vapor  measurerients, 
the  mixing  ratio  is  mostly  given  in  parts  per  million  by  weight  (ppniv.  h 
In  the  present  review,  the  U.  S.  Standard  Atmosphere  [US-62]  has  been 
used  to  convert  between  densities  and  mixing  ratios. 
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Fig.  9-1.  Volume  Mixing  Ratios  for  H2O  Up  To  16  km. 
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The  volume  mixing  ratios  for  water  vapor  at  altitudes  between 
10  and  100  km  are  shown  on  Fig.  9-2. 

The  altitude  profiles  of  mixing  ratios  have  been  derived  from 
infrared  emission  and  absorption  spectra  obtained  by  the  University  of 
Denver  group  on  balloon  flights  |BG-73c,  MG-73a,  GM-73b]  in  February 
and  June  at  midlatitudes.  A minimum  between  15  and  20  km  is  found  and 
this  is  confirmed  by  comparison  of  two  other  optical  measurements 
[Ac-74,  Fa-74a].  The  measurements  of  Farmer  lFa-74a]  in  the  alti- 
tude range  of  14  to  20  km  based  on  the  water  yielded  a con- 

stant mixing  ratio  of  2.4  tO.  3 ppmv  with  no  significant  variation  over 
the  latitude  range  between  30  N and  80  N.  Absorption  measurements  of 
Ackerman  [Ac-74  |,  also  using  the  l^and  of  water,  gave  a mixing  ratio 
of  3.4  ±0.  7 ppmv  in  the  altitude  range  between  20  and  37  km.  The  meas- 
urements of  Ackerman  [Ac-74]  are  also  consistent  with  a slight  and 
broad  maximum  betw’een  20  and  30  km.  The  measurements  of  Goldman 
et  al.  [GM-73b]  also  indicate  a seasonal  variation  in  the  mixing  ratio 
with  summer  values  smaller  than  winter  values.  Other  measurements 
[MG-73a]  by  the  University  of  Denver  group  in  India  showed  a steady  de- 
cline in  mixing  ratio  between  10  and  26  km;  these  measurements  are  not 
shown  on  Fig.  9-2. 

By  means  of  frost-point  temperature  measurements,  a series 
of  midlatitude  measurements  have  been  made  by  Mastenbrook  [Ma-74c, 
Ma-71c|  for  altitudes  up  to  30  km  over  a period  of  ten  years.  These 
median  mixing  ratios  [Ma-74c|  are  shown  in  Fig.  9-2  for  altitudes  be- 
tween 16  and  28  km. 

Data  by  Evans  [Ev-74a],  as  reported  by  Hard  [Ha-74]  but 
without  any  details  of  the  measurements,  are  plotted  on  Fig.  9-2,  having 
been  converted  from  densities  to  mixing  ratio  values. 
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A measurement  of  the  water  vapor  concentration  was  made  at 
28  km  by  use  of  a spin-flip  Raman  laser  technique  [PB-74]  as  con- 
sidered also  in  Section  5.  3.  1;  as  can  be  seen  on  Fig.  9-2,  this  meas- 
urement indicates  a very  dry  atmosphere. 

Emission  measurements  for  water  vapor  using  an  infrared 
radiometer  [HP-74]  have  been  made  in  the  southern  hemisphere  at  34°S 
as  well  as  23°S.  The  October  and  February  measurements  at  34  S are 
shown  in  Fig.  9-2;  the  February  data  are  in  accord  with  other  meas- 
urements, generally,  but  the  October  data  indicate  a wet  layer  around 
20  km.  Emission  measurements  in  the  sub-millimeter  infrared  region 
have  been  made  by  Harries  [Ha-73e]  from  aircraft  at  altitudes  between 
12  and  16  km  over  the  latitude  range  between  50°N  and  40  S.  There  is 
considerably  less  scatter  in  these  data  when  the  volume  mixing  ratio  is 
plotted  against  Ah,  the  difference  in  altitude  between  the  measurement 
height  and  the  tropopause  height,  than  when  the  data  are  plotted  against 
the  measurement  height.  The  data  are  shown  on  Fig.  9-2  as  they  would 
appear  plotted  against  Ah  but  with  a tropopause  height  of  14  km  which 
is  near  the  average  of  the  values,  ranging  between  11  and  17  km,  asso- 
ciated with  the  measurements.  The  data  also  show  no  clear  correlation 
with  latitude.  Harries  et  al.  [HS-73b]  have  also  made  measurements 
of  the  mixing  ratio  for  altitudes  above  15.  2,  24.  3 and  31.  8 km  from  a 
balloon  platform.  The  average  value  was  found  to  be  2.  7 ppm  with  no 
significant  difference  between  the  measurements  at  the  three  altitudes. 
These  measurements  have  not  been  included  in  the  profile  of  Section  9.  4 
since  the  reference  [HS-73b]  was  not  available  in  time. 

Sissenwine  et  al.  [SG-68a]  have  made  17  soundings  over  north- 
ern California  up  to  32  km  using  an  automatic  frost-point  instnament  on 
balloon  flights.  These  data  are  show’n  on  Fig.  9-2  and  exhibit  significantly 
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largei’  mixing  ratio  values  above  20  km  than  any  other  measurements 
shown.  In  the  experiments  [SG-68a]  care  was  taken  to  minimize  water 
contamination,  but  no  quantitative  assessment  of  contamination  appar- 
ently could  be  made.  The  data  profile  has  been  extended  to  80  km  by 
Cadle  lCa-74|  based  on  a typical  temperature  found  in  noctilucent  clouds 
at  80  km.  Sissenwine  et  al.  [Ca-74|  have  considered  convective  clouds 
which  penetrate  the  tropopause  as  a source  of  the  significantly  higher 
water  content  found  above  20  km:  McElroy  (Ca-74  | has  emphasized  the 
source  of  water  resulting  from  oxidation  of  methane  in  the  stratosphere. 
Sissenwine  et  al.  [SG-68a)  also  have  considered  the  passage  of  water 
vapor  into  the  stratosphere  from  lower  altitudes  during  overlapping  of, 
say,  the  tropical  and  polar  tropopauses.  However,  no  definitive  account 
of  the  large  mixing  ratios  found  bV  Sissenwine  et  al.  [SG-68a]  exists. 

From  an  air  sample  collected  over  the  altitude  range  between 
44  and  62  km,  Scholz  et  al.  )SE-70]  have  found  a range  of  mixing  ratio 
values  from  3 to  10  ppmv.  An  analysis  to  distinguish  the  high  altitude 
water  vapor  from  that  desorbed  from  the  sampler  surface  during  anal- 
ysis, had  to  be  made.  The  ranges  of  the  mixing  ratios  and  altitudes  are 
shown  on  Fig.  9-2. 

Attempts  have  been  made  to  determine  water  vapor  densities 
at  altitudes  above  100  km.  One  experiment  [OP-72]  in  which  a mass 
spectrometer  was  used  in  the  altitude  range  between  130  and  210  km, 
was  admittedly  unsuccessful  as  a result  of  water  contamination.  Meas- 
urements of  water  density  in  the  range  from  90  to  180  km  are  shown  by 
Martynkevich  [Ma-72d]  which  are  two  to  four  orders  of  magnitude 
larger  than  values  calculated  by  Martynkevich  [Ma-72d|;  this  excess 
is  attributed  to  hydrogen  and  oxygen  recombination  processes  inside  the 
mass  analyzer.  However,  the  calculations  yielding  the  lowest  water 


density  are  based  on  data  of  Goldberg  and  Blumle  [GB-70]  who  consider 
their  mass  18  data  to  be  entirely  due  to  contamination:  if  this  assign- 
ment is  correct,  the  measured  water  densities  of  Martynkevich  [Ma-72d] 
could  be  as  much  as  four  to  six  orders  of  magnitude  higher  than  calcula- 
tions (see  also  Section  9.3).  Furthermore,  the  mass  18  data  of  Hoffman 
et  al.  [HJ-69a]  which  Martynkevich  also  uses  in  calculating  water  den- 
sity are  attributed  by  the  investigators  |HJ-69a]  to  outgassing  from  the 
rocket  payload  and  this  is  substantiated  by  satellite  measurements  [HJ- 

69a]  not  subject  to  outgassing  and  in  which  the  mass  18  signal  is  attrib- 

18'  2-3 

uted  to  O with  a nuiximum  measured  densitv  of  10  cm  at  about 

700  km. 

9.3  MODEL  CALCULATIONS 

The  results  of  calculations  of  water  density  are  collected  in 
Fig.  9-3.  The  calculations  of  Hunt  [Hu-73a],  Thomas  and  Bowman 
(TB-72|,  Geoi-ge  et  al.  |GZ-72a  ! and  Hunten  and  Strobel  [HS-74a]  are 
compiired.  Between  60  and  80  km,  the  differences  in  the  calculated 
water  densities  are  relatively  small.  Between  80  and  100  km,  the  dif- 
ferences increase  significantly:  these  changes  reflect,  in  a rough  fashion, 
the  relative  values  and  changes  with  altitude  of  the  eddy  diffusion  coeffi- 
cients (K)  used  by  the  calculators.  Thus  the  smallest  values  of  the  eddy 
diffusion  coefficient  in  the  region  from  60  to  about  95  km  were  used  by 
Hunten  and  Strobel  |HS-74a]  whose  calculations  yield  the  smallest  water 
densities  above  80  km.  The  eddy  diffusion  coefficient  profile  used  by 
George  et  al.  [GZ-72a]  has  the  largest  values  of  K at  80  km  but  the 
values  decline  at  higher  altitudes  becoming  le.ss  than  those  used  bv  Hunt 
[Hu-73a]  and  Thomas  and  Bowman  [TB-72]  above  about  92  km  and  less 
than  those  used  by  Hunten  and  Strobel  [HS-74a]  above  95  km.  However, 
the  importance  of  eddy  diffusion  declines  at  these  altitudes  and  becomes 
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Fig.  9-3.  Calculated  Water  Densities  Between  40  and  160  km 


less  impoi-tant  than  molecular  diffusion  at  altitudes  in  the  vicinity  of 
100  km.  In  some  of  the  calculations  [GZ-72a,  TB-72]  the  relative  im- 
portance of  eddy  diffusion  declines  rapidly  at  altitudes  between  95  and 
105  km  whereas  a slower  decline  occurs  in  others  [Hu-73a,  HS-74a]. 

The  relative  values  of  water  densities  established  in  the  four  calcula- 
tions before  reaching  altitudes  slightly  above  100  km  govern  the  differ- 
ences in  the  water  densities  found  at  the  highest  altitudes  of  the  calcu- 
lations. Above  about  100  km,  molecular  diffusion  and  photolysis  of  water 
dominate  the  water  profiles;  as  these  processes  are  similarly  treated 
by  the  calculators,  the  four  profiles  have  only  small  relative  changes  at 
these  highest  altitudes. 

The  results  of  the  four  calculations  [GZ-72a,  HS-74a,  Hu-73a, 
TB-72]  are  also  shown  on  Fig.  9-2  for  comparison  with  the  experimental 
data.  In  addition,  calculations  of  Shimazaki  and  Ogawa  [SO-74c|  for 
altitudes  between  15  and  50  km  and  of  Cadle  (Ca-74  | for  the  extension 
of  the  measurements  of  Sissenwine  et  al.  [SG-68a]  to  80  km  are  plotted 
on  Fig.  9-2.  All  calculations  are  for  midlatitudes. 

9.  4 SELECTION  OF  WATER  PROFILE 

Between  ground  and  16  km,  the  selected  profile  is  given  by  the 
midlatitude  mean  values  assembled  by  Cadle  [Ca-74];  up  to  8 km,  these 
values  are  based  on  the  annual  average  of  radiosonde  data  for  45  N 
latitude  [GS-66bj;  between  10  and  16  km,  these  values  are  based  on 
midlatitude,  a- radiation  soundings  [SG-68a,  SG-68b].  The  data  are 
plotted  in  Fig.  9-4  and  for  altitudes  of  4 km  and  higher;  bars  show  the 
range  of  the  1' o extremes.  The  data  given  by  Cadle  [Ca-74]  have  been 
converted  to  water  densities  by  using  the  U.  S.  Standard  Atmosphere 
[US-62]  total  densities.  See  Table  9-la.  Note  that  the  range  in  density 
of  water  vapor  at  the  surface  of  the  earth  ranges  over  five  orders  of 
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Table  9- la.  Midlatitude  Water  Vapor  Concentrations  in  the  Troposphere 


Average  of  values  over  the  rant;i'  19  to  21  kni. 
Average  of  values  over  the  range  24  to  27  km. 
These  values  rear!  from  HP-74. 


magnitude  and  that  although  the  surface  density  generally  decreases 
with  increase  in  latitude,  there  are  variaticms  due  to  local  sources  and 
sinks  of  water  vapor  [Ca-74  |. 

The  majority  of  experimental  result  >f  Fig.  9-2  lie  between 
15  and  50  km;  between  these  altitudes  th  -ie  i.s  a nominal  range  in  mix- 
ing ratio  between  1 and  30  ppmv.  1 ■ derive  a n-iminal  profile,  the  ex- 
perimental data  in  the  altitude  range  between  15  and  4 5 km  were 
averaged  excluding  the  v^ery  high  ratios  of  Sisscnwine  et  al.  |SG-68a] 
and  the  "wet  layer”  data  of  Hysen  and  Platt  | HP-74]  identified  on 
Fig.  9-2  by  the  filled  circles.  The  data  and  the  averages  are  given  on 
Table  9- lb. 

Between  60  and  80  km  an  average  of  the  calculated  densities 
was  used  in  the  selected  profile.  Above  100  km,  an  average  of  the  cal- 
culated densities  of  Hunten  and  Strobel  [HS-74a|  and  George  et  al. 
[GZ-72a|  was  used  in  the  selected  profile  on  the  basis  of  the  arguments 
of  Hunten  and  Strobel  [HS-74a]  for  lower  values  of  the  eddv  diffusion 
coefficient  (in  the  region  above  80  km).  Thus  the  calculations  of  Hunt 
|Hu-73a]  and  Thomas  and  Bowman  [TB-72|  are  excluded.  .\t  altitudes 
between  4 5 and  60  km  and  between  80  and  100  km,  densities  were  chosen 
by  interpolation.  At  altitudes  above  120  km,  an  exponential  decline  in 
water  density  was  assumed.  Table  9-2  and  Fig.  9-5  show  the  selected 
data.  No  diurnal  variation  was  included  as  the  calculations  ]Hu-73a, 
TB-72]  show  no  significant  variation. 
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